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SUMMARY 


This report describes the results of the data reduction, data analysis, and data archiving of the Pioneer 
Venus (PV) Retarding Potential Analyzer (ORPA) experiment for the period Aug. I, 1988 to Sep 30, 1993. 
The ORPA was one of twelve instruments carried on the PV spacecraft and one of sixteen experimental 
investigations. The Principal Investigator of the ORPA experiment was William C. Knudsen of Knudsen 
Geophysical Research Inc. He has been responsible for representing the interests of the ORPA at PV 
Project Office meetings, directing the data reduction and archiving activities which were performed by 
NASA/Ames, analyzing the reduced ORPA data, and reporting and publishing the scientific results of the 
analysis. In addition, he has provided data, software, and assistance to Prof. Douglas Jones of the Brigham 
Young University, and Dr. Karl Spenner of the German Fraunhofer Institute during approximately the last 
two years of the contract period. Prof. Jones and Dr. Spenner are analyzing ORPA solar wind electron data 
and entry ionospheric data, respectively. 

Periapsis altitude of the PV spacecraft descended slowly from approximately 1900 km in Aug. 1988 to 
approximately 130 km in Oct. 1992 at which time communication with the spacecraft ceased. The 
spacecraft was above the ionosphere for most of the contract period, dipping into the ionosphere in the 
nightside hemisphere primarily during the last three months prior to loss of the spacecraft. Because of this 
sampling restriction, the ORPA was commanded into its suprathermal electron mode during most of the 
contract period, a mode in which it measured solar wind electron data, and was commanded into its 
ionospheric measuring modes only when the spacecraft began dipping into the ionosphere toward the end 
of the mission. 

The quantity of data recorded during the last, approximately four years of the PV spacecraft lifetime 
covered by this contract was reduced considerably from that obtained previously because of the reduction 
in power available for operating the science instruments as a result of spacecraft solar cell deterioration, 
loss of half of the spacecraft solid state memory for storing science data during occultations, and partial 
failure of the ORPA during the last two years of spacecraft life. 

Despite the reduction in available data, especially ionospheric data, important scientific discoveries 
resulted from this contract. The most significant was the discovery of a strong solar cycle change in the 
size of the dayside ionosphere and the resulting ’shutoff of flow of dayside ions into the nightside 
hemisphere. The large topside 0+ F2 ionospheric layer observed in both the day- and nightside 
hemispheres during the first three years of the Pioneer Venus mission is essentially absent during the solar 
minimum activity period. 


INTRODUCTION 


The Pioneer Venus retarding potential analyzer instrument (ORPA) is one of eleven instruments that was 
flown on the Pioneer Venus spacecraft, a spacecraft that was launched May 20, 1978 by the United States 
National Aeronautics and Space Administration (NASA), was inserted into orbit about Venus on December 
4, 1978, and ceased transmitting data to earth following the periapsis pass of orbit 5056 The last data 
recovered was from a periapsis period from orbit 5055 on Oct. 7, 1992. The purpose of the Pioneer Venus 
(PV) Mission was to significantly increase our knowledge of the plasma environment, atmosphere, surface, 
and interior of planet Venus. 

The ORPA was selected for the PV mission to measure several important ionospheric plasma quantities as 
the spacecraft swept through the ionosphere near the orbit periapsis. Quantities measured uniquely by the 
ORPA include ion temperature, ion bulk velocity (vector), and suprathermal electron temperature. 
Quantities measured in common with other instruments include total ion (plasma) density, major ion 
concentrations, electron temperature, and spacecraft potential. The ORPA has also proven capable of 
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measuring the solar wind electron temperature and density, a capability which has proven to be of great 
value in defining plasma properties in the boundary layer between the solar wind and ionosphere. 

The PV mission data acquisition, reduction, analysis and dissemination activities since December 4, 1978 
have been separated into four programmatic time periods, Prime Mission, May 1978 - Sept. 1979; 
Extended Mission I, Oct. 1979 - Mar. 1981; Extended Mission II, Mar. 1981 - Sept. 1986; and Extended 
Mission III, Sept. 1986 - Sept. 1993 The data reduction and analysis activities reported herein began Aug. 
1, 1988 and extended through Sept. 1993 and covered most of the Extended Mission III period. For 
scientific purposes, the mission has been divided into three mission phases. Phase I began at the time of 
insertion of the PV spacecraft into orbit about Venus and ended to approximately Dec. 1980 at which time 
spacecraft fuel ceased to be used to maintain periapsis altitude at low altitude (Fig. 1). Phase II extended 
from approximately Dec. 1980 to April 1991. Phase III extended from April 1991 to termination of 
communication with the spacecraft which occurred on Oct. 8, 1992. The data acquisition period covered 
by this final report is from Aug. 1 1988 through Oct. 7 1992. 

Solar conditions changed during this data acquisition period from solar cycle minimum at the beginning 
(1988) to solar cycle maximum, 1989-1990, and back down toward solar cycle minimum at the end. 
Periapsis altitude decreased from its maximum value of approximately 2300 km at the beginning of the 
contractual period to altitudes as small as approximately 130 km during the final low altitude maintenance 
period. This combination of conditions prevented in situ measurements of the ionosphere during solar 
minimum conditions existing during phase II but provided a short period of in situ measurements at low 
altitude in the nightside ionosphere during solar conditions approaching that of solar minimum during the 
final entry period of the mission. ° 

Despite the paucity of in situ measurements of the ionosphere during solar minimum conditions, the 
measurements that were made provided the clues and, later, the evidence for the "breathing" of the Venus 
ionosphere that occurs through a solar cycle. This breathing was not anticipated at the end of the prime 
mission and is probably the most significant contribution to the overall mission that was made during the 
contractual period covered by this final report. 

The ORPA data analysis and reporting activities for the subject contract have been performed by William 
C. Knudsen (WCK). The data reduction has been performed at NASA/Ames Research Center under the 
supervision of WCK with Joe Goosby, a contract employee with Sterling Software performing 
programming and computer operations. 


MISSION OPERATIONS 


The mission operations during this contract affected ORPA mode commanding and data recovery in the 
following way: The declining power output of the spacecraft solar cells dictated that the ORPA be turned 
on during a limited portion of each orbit. Hence, it was necessary to turn the ORPA off and on at least 
once each orbit. During the early years of the mission, the ORPA was left on at all times, and commanding 
it consisted of sending mode reconfiguration commands only. The spacecraft data storage unit lost half its 
data storage capacity which cut in half the data that could be stored during short occultations (periods when 
the spacecraft was shadow from the Earth by Venus). The orbit periapsis altitude was above the 
ionosphere except for the last year of the mission. The consequences of these conditions were that the 
ORPA was commanded into its suprathermal electron mode most of the time except for the last year. 
During the last year, It was commanded into its thermal electron and ion modes during short periods 
around periapsis. Because of the limited power output of the solar cells, the loss of half the data storage 
unit, competition among principle investigators for data formats and bit rates, and a partial failure of the 

ORPA, ionospheric data obtained by the ORPA during the final eclipse season of the mission was severely 
limited. 
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ORPA HEALTH 


The ORPA operated within its specifications until June 1990 at which time it began to respond to 
commands erratically. The data stream from the instrument was observed on some occasions to be 
anomalous after turn-on and mode configuration by ground command. At this phase of the mission, the 
ORPA, as well as other instruments, was turned off during portions of the orbit to conserve power. After 
experimentation over a period of approximately one year, I determined, with the cooperation of the Pioneer 
operations center, that the ORPA was accepting ground commands correctly approximately 50% of the 
time. The remainder of the times, the ground commands were implemented incorrectly, and the bit stream 
from the instrument was anomalous. A procedure was implemented with the operations center such that 1 
the bit stream was observed to be anomalous after ORPA tum-on, the ORPA would be turned off and then 
re-commanded. Because of the two-way time delay between sending a command and receiving the 
instrument response, this procedure was possible only once per orbit and, then, not on all orbits. With this 
procedure, we achieved approximately 70% data recovery of the total possible recovery during the last year 
of spacecraft operation. 


PROJECT OFFICE SUPPORT 


WCK attended all Science Study Group meetings conducted during this contract period and presented the 
results of his analyses therein. He also attended all of the Orbiter Mission Operations Planning (OMOP) 
meetings and contributed to the deliberations as they affected the ORPA. 


DATA REDUCTION 


The reduction of the ORPA bit stream to engineering numbers and, thence, to geophysical quantities such 
as ion temperature and density, the computer manipulation and plotting of the geophysical quantities to aid 
analysis, and the submission of the geophysical quantities in appropriate format to the National Space 
Science Data Center (NSSDC) was performed at NASA/Ames Research Center under the direction of a 
NASA appointed PV ORPA data manager. The PV ORPA data manager was Robert C. Whitten, initially, 
and, thereafter, John Mihalov. The programming and computer operations required to accomplish the 
reduction were performed primarily by Mr. Joseph Goosby, a contract employee of Sterling Software. 

The programming and data reduction activities performed by Mr. Goosby has been disappointing. I have 
been unable to accomplish all the analysis activities that I could have accomplished had the programming 
and data reduction been performed more effectively. I have had to restrict my requests for data analysis 
because of delays and lack of accomplishing the programming and routine data reduction. 

All of the ORPA data from orbit 1 to orbit 5055 have been re-reduced during this contract period with 
revised reduction software programs that were written by WCK to improve the accuracy of the derived 
geophysical quantities. With the exception of the thermal electron mode, derivation of geophysical 
quantities from the I-V curves returned by the ORPA requires that the 1-V curves be fit with a theoretical 
expression that predicts the ion or electron current that should be observed as a function of retarding 
potential. The theoretical expression depends on instrument parameters as well as the ambient plasma 
geophysical parameters. WCK spent a substantial fraction of his efforts during this contract modifying and 
improving the theoretical expressions for the ion current when the ORPA is operating in its ion mode and 
for the electron current when it is operating in the suprathermal electron mode. 
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NSSDC DATA SUBMISSIONS 


The first product that has been submitted to the National Space Science Data Center (via the Planetary Data 
System (PDS) Node at the University of California) is called the RDR files and contains the ORPA current- 
voltage (I-V) data for all retarding potential sweeps for all orbits. The I-V data are in engineering 
quantities, Amperes and Volts, and are associated with a status word which describes the instrument mode 
in which the data were taken and a set of ephemeris quantities which describe the time of measurement and 
the orientation and location of the instrument and spacecraft. The RDR files will be essential for verifying 
the accuracy of the least-squares fits from which the geophysical quantities to be submitted as the Tables 
files and UADS files were derived. PDS has transcribed the RDR files from 9 track magnetic tapes to CD 
ROMS for convenient use by the scientific community. The RDR files for the final entry period of the PV 
mission are presently being used by Dr. Karl Spenner of the German Fraunhofer institute in Freiburg to 
derive the geophysical quantities for the entry period. Dr. Spenner was a Co-investigator on the ORPA 
experiment during the prime PV mission period, and was responsible for supplying the sensor head for the 
ORPA. 

Two additional products which are to be submitted to NSSDC are called the Tables File and the UADS 
File. The Tables File contains all the geophysical quantities derived from the ORPA for the entire PV 
mission. It is complete and archived at NASA/Ames Research Center at the present time. The UADS File 
will contain ORPA geophysical quantities for +/- one half hour about periapsis for each orbit and will 
conform to the format specified by the PV Science Study Group for this File. The UADS file contains a 
subset of the geophysical quantities stored in the Tables file. The PV ORPA data manager has been 
instructed to submit the ORPA UADS file consistent with the "letter to distribution" by Lawrence Colin 
dated Sept. 7, 1983, the report referred to therein and titled "Pioneer Venus Mission, Instructions for Data 
Submissions to the National Space Science Data Center, A Committee Report" by Roger Craig, NASA 
Ames Research Center dated Aug. 11 1983, and the document "Pioneer Venus Orbiter Retarding Potential 
Analyzer NSSDC Submission Documentation" by William C. Knudsen Dated Nov. 4, 1993. The letter, 
report and document are reproduced in Appendix C. Submission of the ORPA reduced data in the UADS 
form satisfies the contract requirement for data submission to NSSDC. 


ANALYSIS RESULTS 


Improvement of Reduction Algorithms 

A substantial fraction of WCK's efforts during this contract has been expended on improving the least- 
squares algorithms used to obtain geophysical quantities from the suprathermal electron and ion modes of 
the ORPA. The algorithms used prior to the re-reduction of all orbits, which took place in the Spring of 
1993, were based on the assumption that the grid walls of the ORPA were infinite in radius and the grid 
thicknesses were zero. He believed that these assumptions were most likely to be affecting the density and, 
to a smaller extent, the temperature of the suprathermal electrons derived from the I-V curves. In an early 
publication containing ORPA solar wind electron densities and orbiter plasma analyzer (OPA) solar wind 
ion densities, it was evident that the ORPA densities were consistently a factor of 2 below the OPA ion 
densities. To try and improve the agreement of the densities from the two instruments, he undertook a 
numerical analysis of the reduction of the electron current to an RPA having the geometry and grid 
thicknesses of the ORPA caused by the existence of a finite wall radius and finite thickness of grids. The 
results of this analysis have been published [Knudsen, 1992a] and incorporated into the ORPA data 
reduction program. The new algorithm is now providing solar wind densities that are approximately a 
factor of two larger than those obtained without the correction for finite wall radius and grid thickness. 
Professor Douglas Jones, of the Brigham Young University, is currently analyzing solar wind electron 
density and temperature data from the ORPA and will publish a detailed comparison of densities derived 
from the ORPA and OPA. 
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Solar Cycle Changeiia i hs. . Ve nus Ionosphere 


The most significant achievement of this contract effort is the discovery by WCK of a rather dramatic solar 
cycle variation in the size of the Venus ionosphere. The Venus ionosphere, as initially measured early in 
the PV mission, typically extended from below 150 km altitude to 300 km altitude at low solar zenith 
angles (SZA) and to 900 km at 90 deg SZA. In the nightside hemisphere, an upper ionospheric layer of 
atomic ions existed which extended over 900 Km in altitude and had a density of the order of lx 10+3 
cm-3. This was the picture generally held of the Venus ionosphere prior to his discoveries reported in 
publications [Knudsen, et al. 1986; Knudsen et al. 1987; Knudsen, 1988; and Knudsen, 1992b]. A 
summary of the Venus ionosphere as understood up through 1991 is presented in the publication "The 
Venus Ionosphere from In-situ Measurements" [Knudsen, 1992b] and is included in this report as 
Appendix A. The solar cycle variation as described therein is now generally accepted and has been 
verified by the reentry data obtained in the last 50 orbits of the PV mission and which are currently being 
reported in presentations and publications by several of the PV principal investigators. 

Additional Ionospheric Results 

Two additional papers on ionospheric topics, "Role of Hot Oxygen in Venusian Ionospheric Ion Energetics 
and Supersonic Antisunward Flow" [Knudsen, 1990a] and "The Venus Transterminator Ion Flux at Solar 
Maximum" [Knudsen and Miller, 1992], were also published during this contract. The first paper drew 
attention to the role of hot oxygen in heating the Venus ionosphere and in affecting the transterminator 
flow of dayside 0+ ions. The second paper presented a more detailed and accurate measurement of the 
average transterminator flux of 0+ ions during solar maximum activity than had been presented in an early 
paper. In the second paper, evidence was again presented in favor of the belief that the major source of 
ionization of the nightside ionosphere during solar maximum activity is transterminator 0+ transport. 

Supr^emd.El ect ron M easurements 

One of the important goals of this contract was to obtain as much information from the suprathermal 
electron mode measurements as possible since the PV orbit periapsis was above the sensible ionosphere 
most of the contract period. Valuable suprathermal electron data could be obtained in the mantle, in the 
solar wind in the dayside hemisphere, and in the ionotail (wake of the planet). 

Possible mass loading of the solar wind up stream of the Venus bow shock has been an interesting subject 
because of the relatively close proximity of the Venus bow shock (approximately.2000 km above the 
Venus surface in the sub solar region). Brace et al. [3985] and Brace [1987] observed increased electron 
current to the electron temperature probe (ETP) upstream of the Venus bow shock in the sub solar region 
and referred to this region as the precursor region. They suggested that mass loading of the solar wind 
upstream of the bow shock through charge exchange and photo ionization of the neutral, hot atomic 
oxygen halo was responsible for the increased current. Knudsen [1985] presented observations by the 
ORPA of a smooth increase in solar wind electron density and temperature in the same region and assumed 
the cause to be mass loading of the solar wind. In a subsequent paper titled "The Venus Precursor: An 
Environmental Effect on the Pioneer Venus Spacecraft", Knudsen et al. [1989] concluded that both the 
increase in current observed by the ETP and the increase in density of a relatively cold (about 1 eV in 
temperature) electron population measured by the ORPA were caused by an increase in the density of 
secondary electrons trapped in the potential well surrounding the PV spacecraft. The increase in Solar 
wind electron density and temperature which had been reported by Knudsen [1985] was shown to be 
caused by the grazing passage of the PV spacecraft into and through the bow shock transition region. 
Knudsen et al. [1989] concluded, further, that there was no evidence in their data for a mass loading of the 
upstream solar wind density that exceeded a few percent. 

The source of the suprathermal electrons observed in the Venus umbra with an energy distribution 
characterized by a temperature of about 14 eV has been an intriguing question ever since the electrons were 
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first measured by Gringauz et al. [1979] with a retarding potential analyzer carried on the Venera 9 
spacecraft in 1976. The same population has been measure by the PV ORPA, and the source of the 
electrons has remained unknown. However, measurements of the suprathermal electron temperature and 
density as the PV spacecraft travels from the magnetosheath into the planet umbra strongly suggest that the 
source of the electrons is the solar wind. The temperature of the suprathermal electron distribution 
decreases smoothly from that of the shocked solar wind flowing in the magnetosheath to that of the 
suprathermal electrons co-located with the ionospheric plasma in the umbra. Frequently, but not always, 
the density decreases in like manner. This variation of the suprathermal electrons with location strongly 
suggests that the suprathermal electrons are solar wind electrons that have infiltrated the nightside 
ionosphere and ionotail. This suggestion and the measurements supporting it have been presented at an 
annual meeting of the American Geophysical Union [Knudsen, 1990b]. A paper containing these results 
will be published in due time. 

The results of a study of the suprathermal electron flux in ionospheric holes are presented in Appendix B. 
The results, titled "Suprathermal Electron Flux in Venus Ionospheric Holes" was submitted for publication 
but was rejected, because only two examples were presented. From these two examples and additional 
examples not presented, WCK concluded that the suprathermal electron flux observed in ionospheric holes 
is similar to that observed outside the holes and that the holes are not formed by the existence of a 
significant electric field parallel to the field as has been by suggested Grebowsky and Curtis [1981]. The 
cause or causes of ionospheric holes remains a significant unsolved problem, and WCK intends to pursue 
this problem by examining additional examples. 

WCK contributed data and analysis to two studies undertaken by Dr. Devrie Intriligator. One study 
resulted in the publication of a paper "Evidence for Unusually high densities of Plasma in the Venus 
Ionosheath" [Intiligator et al., 1991], This study attributed the high density of hot plasma to ionospheric 
plasma that had been drawn into and heated in the ionosheath. The second study, "Evidence for Ion 
Transport on the Venus Ionotail" is an effort to report on the first believable low energy ion fluxes 
measured by the Plasma analyzer and relate these measurements and related measurements by the other 
in situ instruments to a narrow plasma feature observed on orbit 4574. This paper has not been accepted 
for publication at the present time. 


PV SPACECRAFT ENTRY RESULTS 


Descent of the PV orbit periapsis altitude into the ionosphere during the last few weeks prior to spacecraft 
bumup was a particularly important event of the PV mission. The periapsis passes occurred in the 
nightside ionosphere during a period of time when the Sun's activity was approaching a minimum, and the 
periapsis altitude was permitted to get smaller than it was during the early years of the mission since the 
risk of losing the spacecraft was of little consequence. The significance of the entry occurring at a time of 
near-minimum solar activity is that in situ measurements could confirm predictions that the transport of 
dayside 0+ ions across the terminator would be reduced a factor of ten during solar minimum. WCK and 
his co-authors had predicted in 1986 and 1987 that the nightside topside 0+ layer density would be 
reduced a factor of ten as a result of the reduction in transterminator 0+ flow [Knudsen et al., 1986; 1987]. 
Preliminary analysis of a few ORPA orbits and reported results from the orbiter ion mass spectrometer and 
electron temperature probe instruments for the entry period confirm the accuracy of the prediction. 

Unfortunately, the amount of ORPA data obtained during this period was severely limited by the reduced 
interval of time that the spacecraft could store data, the limited bit rate, and the partial failure of the ORPA. 
Nevertheless, significant results will be forthcoming from the ORPA entry data. Analysis of the ORPA 
entry data may not be accomplished under this contract because of slow progress in processing of ORPA 
data at Ames Research Center. If no results are forthcoming from this contract, results will be forthcoming 
from an analysis effort by Dr. Karl Spenner sponsored by the German government. The effort is a three 
year effort and is in support of analysis of ORPA data obtained during the entry period. Dr. Spenner has 
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been supplied our data reductions programs and the entry ORPA data. He is currently reducing the data by 
interactively performing least-squares fit of the theoretical ion current algorithm, mentioned above, to the 
entry ORPA I-V data. Interactive fitting of the algorithm to the data is essential because of the meager 
amount of ORPA data and the low level of the currents. 


PUBLICATIONS AND PRESENTATIONS 


Publications 

The scientific papers written and completed during this contractual period are: 

Phillips, J. L., J. G. Luhmann, W. C. Knudsen, and L. H. Brace, Asymmetries in the location of the Venus 
ionopause, J, Geophys. Res., 93, 3927-3941, 1988 

Knudsen, W. C., Solar cycle changes in the morphology of the Venus ionosphere, J. Geophys. Res., 93, 
8756-8762, 1988. 

Knudsen, W. C., J. G, Luhmann, C. T. Russell, F. L. Scarf, The Venus precursor: An environmental effect 
on the Pioneer Venus spacecraft, J. Geophys. Res.,94, 197-207,1989. 

Knudsen, W. C., Role of hot oxygen in Venusian ionospheric ion energetics and supersonic antisunward 
flow, J. Geophys. Res., 95,1097-1101,1990. 

Singhal, R. P., R. C. Whitten, and W. C. Knudsen, Comparative dynamics of the ionospheres of Venus and 
Mars at large solar zenith angles, Indian J. Rad. and Space Sci.,9, 88-102, 1990. 

Intriligator, D. S., L. Brace, S. H. Brecht, W. C. Knudsen, F. L, Scarf, R. J. Strangeway, and H. Taylor, 
Evidence for unusually high densities of plasma in the Venusian ionosheath, Geophys. Res. Let., 18, 61-64, 
1991. 

Knudsen, W. C., Finite grid radius and thickness effects on RPA- measured suprathermal electron density 
and temperature, J. Geophys. Res., 97, 13,767, 1992a. 

Knudsen, W. C., The Venus Ionosphere from In-situ Measurements, in Venus and Mars: Atmospheres, 
Ionospheres, and Solar Wind Interactions, J. G. Luhmann, M. Tatrallyay, and Robert O. Pepin, Eds., pp. 
237 - 263, AGU Geophysical Monograph 66, The American Geophysical Union, Wash. DC., 1992b. 

Knudsen, W. C. and K. L. Miller, The Venus transterminator ion flux at solar maximum, J. Geophys. Res., 
97, 17,165, 1992. 

Knudsen, W. C., Venus' internal magnetic field and its reconnection with the IMF, J. Geophys. Res., 
submitted 1993. 

Presentations 

The presentations of scientific results derived from the data reduction and analysis of ORPA and other PV 
instruments are as follows: 

Knudsen, W. C., Current picture of the Venus ionosphere from Pioneer-Venus and major unknowns, 
Twenty-seventh Plenary Meeting of COSPAR, ABSTRACTS, p237, Espoo, Finland, July, 1988. 
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Knudsen, W. C., Role of hot oxygen in Venusian ionospheric ion energetics and antisunward flow, EOS\ 
70, 1177, 1989 

Knudsen, W. C., The Venus Ionosphere from In Situ Measurements, AGU Chapman Conference on 
Comparative Study of Venus and Mars, Balatonfuered, Hungary, June 3-9, 1990. 

Knudsen, W. C., Origin of Suprathermal Electrons in the Venus Ionosphere, EOS, 71 , 1431, 1990. 

Knudsen, W. C., Pioneer Venus Mission and Venus Plasma Environment, presented at the April 23 
American Institute of Aeronautics and Astronautics, San Francisco Section, April 23, 1992a. 

Knudsen, W. C., Venus internal magnetic field and its interaction with the interplanetary magnetic field, p 
60, International Colloquium on Venus, Aug. 10-12, 1992b, Pasadena, CA. 

Knudsen, W. C., Venus 1 Internal Magnetic Field and its Interaction with the Interplanetary Magnetic field, 
EOS, Trans., Amer. Geophys. Union, 74 , 111, 1993. 


CONCLUSIONS 


The PV mission has been an outstanding success, and the results from the ORPA have contributed greatly 
to that success despite the fact that it had a small -probably the smallest- share of the spacecraft bit stream. 
Significant experimental contributions made by the ORPA include planet-wide measurement of the 
ionospheric electron and ion temperature fields, ion bulk velocity field, total and major ion densities, 
ionopause altitude, and photoelectron energy distribution. The ORPA has also contributed suprathermal 
electron density and energy distribution measurements in the region immediately adjacent to the 
ionosphere but dominated by the solar wind plasma. 

Major scientific discoveries resulting from ORPA measurements and published by ORPA team members 
include: supersonic bulk flow of the ionospheric plasma across the Venus terminator, acceleration of the 
supersonic flow by the thermal pressure gradient of the ionosphere across the terminator, maintenance of 
the nightside ionosphere during solar activity maximum conditions primarily by the supersonic flow of the 
dayside ionospheric plasma into the nightside hemisphere, existence of a transition region named the 
mantle above the ionopause, and compression of the ionopause to low, median altitude by the solar wind 
dynamic pressure during solar activity minimum conditions with consequent large reduction in the density 
of the nightside 0+ layer. 
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APPENDIX A 

THE VENUS IONOSPHERE FROM IN SITU MEASUREMENTS 

BY 

WILLIAM C. KNUDSEN 
Reproduced from: 

Venus and Mars: Atmospheres, Ionospheres, and Solar Wind Interactions, J. G. Luhmann, M. Tatrallyay, 
and Robert 0. Pepin, Eds., pp 237-263, AGU Geophysical Monograph 66, The American Geophysical 

Union, Washington D.C. 1992. 
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THE VENUS IONOSPHERE FROM IN SITU MEASUREMENTS 


zt 


I 


William C. Knudsen 

Knudsen Geophysical Research, Inc., Monte Sereno, California 

A current picture of the Venus ionosphere is presented with emphasis placed on those features 
revealed by the Pioneer Venus orbiter (PVO) in situ ionospheric particle measuring instruments. The 
PVO mission has extended over a complete cycle of solar activity and has revealed a striking change 
in the typical ionospheric morphology from solar maximum to solar minimum. However, the 
periapsis altitude history has been such that the ionosphere was sampled by the PVO in situ 
experiments only during periods of solar maximum. At solar maximum the typical ionosphere, in 
both the day and nightside hemispheres, consists of a lower Fj-type layer composed primarily of O^ 
ions overlain by a substantial F 2 -type layer composed primarily of 0 + ions. The dayside 0 + layer is 
terminated above in a narrow layer called the ionopause across which the plasma origin changes from 
that of the planet Venus to that of the sun. The median altitude of the dayside ionopause at solar 
maximum increases from approximately 300 km to 900 km as the solar zenith angle (SZA) increases 
from 0 to 90° and is determined by pressure balance between the ionospheric particle pressure and 
the normal component of the solar wind dynamic pressure for low to moderate values of the dynamic 
pressure. The ionopause descends to a lower altitude limit of approximately 220 km for high solar 
wind dynamic pressure. The principal source of ionization in the dayside ionosphere during both 
solar maximum and solar minimum is solar extreme ultraviolet (EUV) radiation. The principal source 
of ionization in the nightside ionosphere during solar maximum is a downward flux of 0 + ions 
originating in the day side ionosphere and transported across the nightside hemisphere by a supersonic 
ion wind. The Venus electron temperature is hotter than that of the earth’s ionosphere, exhibits little 
variation with SZA, and is believed to require an external source of heat such as heat flow through the 
ionopause for its explanation. The ion temperature is also hotter than that of the earth's ionosphere 
and is also believed to require an external heat source such as heat flow through the ionopause. At 
higher altitudes in the central nighside hemisphere, a substantial increase in ion temperature occurs 
and is believed to result from conversion of directed bulk ion kinetic energy into thermal energy. 
Narrow regions of low density plasma and enhanced magnetic field directed parallel to the sun-Venus 
axis, called holes, are frequently observed in the nightside ionosphere. The dayside ionosphere 
contains fluxes of suprathermal electrons with energy spectra consistent with the fluxes originating 
from EUV photo production. Temporally and spatially variable suprathermal electron fluxes are 
observed in the nightside hemisphere during both solar maximum and solar minimum with a median 
energy spectrum similar to a 14 eV thermal distribution. These suprathermal electron fluxes are the 
principal ionization source of the nightside ionosphere during solar minimum and are responsible for 
the observed nightside aurora. Although the ionosphere was not sampled by in situ experiments 
during solar minimum, radio occultation measurements and inferences drawn from solar maximum 
in situ measurements reveal a solar minimum ionosphere that consists primarily of an F,-type 0 2 + 
layer in both the day- and nightside hemispheres. The ionopause altitude is approximately 250 km or 
less at all SZA, and the F 2 -type 0 + layer is missing. 


Introduction 

The Venus ionosphere has been sampled for over 12 years by 
several particle measuring instruments mounted on the NASA 
Pioneer-Venus orbiter (PVO) spacecraft [Colin, 1980] and, as a 

Venus and Mars: Atmospheres, Ionospheres, 
and Solar Wind Interactions 
Geophysical Monograph 66 
Copyright 1992 American Geophysical Union 


consequence, is the most extensively explored ionosphere in the 
Solar System next to that of the earth. The intent in this paper is to 
present a current picture of the Venus ionosphere as revealed by 
those instruments. The presentation would be incomplete without 
some reference to important features revealed by the PVO 
magnetometer and radio occultation experiments and, also, radio 
occultation and in situ particle experiments carried on Russian 
Venera missions prior to insertion of the Pioneer Venus Orbiter 
spacecraft into orbit. Previous reviews of the Venus ionosphere 
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covering some of the same material include Brace et al. [1983a], 
Bauer et al. [1985], Luhmann [1986], Brace and Kliore [1991], 
Miller and Whitten [1991], and Luhmann and Cravens [1991]. 

The PVO ionospheric particle instruments include the electron 
temperature probe (OETP), the ion mass spectrometer (OIMS) and 
the retarding potential analyzer (ORPA). The OETP consists of two 
cylindrical Langmuir-type probes, each located at the end of a boom 
approximately one meter long. The instrument measures electron 
temperature, electron density, and total ion density [Krehbiel et al., 
1980]. The OIMS is a radio frequency Bennett ion mass 
spectrometer capable of resolving ions differing in mass by one 
atomic mass unit (amu). The mass range covered is 1 to 56 amu. 
The OIMS also measures the component of the ion drift velocity 
parallel to the instrument axis which is parallel to the spacecraft 
spin axis. When the kinetic energy of the ions relative to the 
spacecraft becomes large, 20-90 eV, the ions are considered 
superthermal, and the density, masses, and kinetic energy of the 
ions are inferred from a harmonic peak analysis method [Taylor et 
al., 1980a, 1980b]. The ORPA is composed of several parallel. 


plane grids followed by a collector. By applying appropriate 
control voltages on the grids, electrons or ions can be selected for 
energy analysis. The ORPA measures electron temperature, ion 
temperature, major ion composition, ion vector velocity, and the 
suprathermal electron energy distribution [Knudsen et al., 1980c]. 

Overview 

An overview of Venus's ionospheric morphology at solar 
maximum is illustrated schematically in Figure 1 together with 
some of the processes that have been observed. As a First 
approximation, the ionosphere may be considered symmetric about 
the sun-Venus axis (the x axis, positive toward the sun in Figure 1). 
A bow shock is formed upstream of Venus in the supersonically 
flowing solar wind and the flow is diverted around the planet. The 
idealized surface separating the planetary ionosphere and the solar 
wind is called the ionopause. The region between the bow shock 
and the ionopause is referred to as the ionosheath or, alternatively, 
the magnetosheath. Just above the ionopause is a region called the 
mantle, a region in which the properties of the ionosheath are 



Solar Maximum 

Fig. 1. Schematic representation of the Venus plasma environment during solar maximum. TTic identified plasma subregions are 
drawn approximately to scale, and several identified processes affecting the the plasmas are shown in their approximate location. 
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altered by the presence of neutral atoms extending above the 
ionopause. The ionosphere itself extends from the ionopause down 
to an altitude of approximately 130 km. The main or peak 
ionospheric layer is predominantly C> 2 + and lies between altitudes 
of approximately 130 and 180 km in the dayside hemisphere and 
between 130 and 160 km altitude in the nightside hemisphere. 
Above the 0 2 + layer, 0 + ions are the dominant constituent in the 
dayside hemisphere. The 0 + and H + ion densities are comparable 
in the nightside hemisphere. The median locations of the bow 
shock and ionopause in Figure 1 are drawn approximately to scale. 

Some of the prominent features and processes characterizing the 
ionosphere at solar maximum include waves on the ionopause, 
elongated Filaments of plasma attached to the ionosphere proper, 
clouds of plasma above the ionopause, the acceleration of dayside 
ionospheric plasma across the terminator into the nightside 
hemisphere at supersonic speeds, the generation of plasma waves 
nightward of the terminator at low altitude, the formation of 
relatively narrow magnetic tubes of low density plasma, called 
holes, oriented approximately parallel to x axis, and the 
precipitation of approximately 14 eV suprathermal electrons into 
the nightside neutral atmosphere, producing ionization and auroral 
emissions (see the papers by Fox and Paxton, this volume). 


The present picture of the solar minimum Venus ionosphere is 
illustrated in Figure 2. Figure 3 shows that the PVO spacecraft 
periapsis altitude has been above 1800 km during the years of solar 
minimum. For this reason, we have no in situ measurements of the 
main ionospheric layer during this phase of the solar cycle so that 
the picture in Figure 2 is based largely on extrapolation and remote 
sensing. The dimensions of the solar minimum main layer has been 
revealed primarily by radio occupation electron density profiles 
recorded during solar minimum (see the paper by Kliore , this 
volume). However, some inferences concerning the composition 
and dimensions of the ionosphere at solar minimum have been 
made from in situ measurements made at solar maximum during 
periods of extreme solar wind dynamic pressure. 

The solar minimum dayside layer is expected to be composed 
primarily of 0 2 + ions and is confined between approximately 130 
and 300 km altitude. The nightside layer is also believed to be 
composed primarily of C> 2 + ions and is confined between 
approximately 130 and 200 km altitude. The in situ measurements 
made above 1800 km reveal that the wake consists of what are 
presumed to be elongated filaments of ionospheric plasma with a 
density of a few times I0 2 cm' 3 . 0 + is a major ion constituent, and 
H + , although not directly measured, is also presumed to be a 
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Fig. 2. Schematic representation of the Venus plasma environment during solar minimum. The PVO orbiter periapsis altitude was 
above 1900 km, and in situ observations of the ionosphere were limited to the relatively high altitude wake region. 
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significant constituent. Between the filaments the density is of the 
order of 10 cm a value comparable to that of the ionosheath in the 
wake. The magnetic field within the troughs is parallel to the x axis 
and is stronger than it is within the filaments. A mantle region 
should be present, but its existence has not been determined by 
measurement. Although Figures I and 2 illustrate what is believed 
to be the typical morphology of the Venus ionosphere at solar 
maximum and solar minimum, respectively, the solar maximum 
ionospheric morphology is believed to look like that of Figure 2 
during periods of extreme solar wind dynamic pressure, and the 
solar minimum ionospheric morphology can presumably look like 
that of Figure 1 during periods of very low solar wind dynamic 
pressure. 

The Solar Maximum Ionosphere 

Ionopause 

One of the Venus ionospheric structures that has no counterpart 
in the ionosphere of the earth is the relatively narrow transition 
region between the solar wind plasma and the ionospheric plasma 
called the ionopause. Figure 4 illustrates the variation in density of 
0 + and 0 2 + ions measured by the PVO ion mass spectrometer as the 
spacecraft passed through the ionopause at a time when the solar 
wind dynamic pressure was low (1.8 nanoPascals (nPa), where 
InPa = IxlO* 8 dynes cm* 2 ) relative to (he median effective 
ionospheric peak particle pressure of 6 nPa, and at a time when the 
solar wind dynamic pressure was high (8.5 nPa) relative to the peak 
ionospheric particle pressure [Mahajan and Mayr, 1989]. Both of 
these orbits were recorded under solar maximum conditions in the 


subsolar region. At low solar wind dynamic pressure, a thick F 2 
type O* layer overlies the layer, and the transition between the 
solar wind plasma and the 0 + layer is characterized by a drop in 0 + 
density by two to three orders of magnitude over an altitude interval 
of a few tens of kilometers. The density scale height within the 
ionopause is small in comparison with the approximately diffusive 
equilibrium 0 + scale height just below the ionopause. In this type 
of transition, the total pressure just above the ionopause is 
principally magnetic field pressure and just below the ionopause is 
principally thermal particle pressure. This thin type of ionopause 
was the most frequent type during the years of solar maximum in 
1979-1981, a period during which the PVO periapsis altitude was 
below 500 km altitude. For large dynamic pressure, the O* layer is 
eroded down upon the 0 2 + layer, and the density scale height within 
the transition layer is larger than that in the low dynamic pressure 
case. The ionospheric plasma density scale height is approximately 
equal to the neutral atomic oxygen scale height [Elphic et a!., 1981; 
Shinagawa et a!., 1987; Mahajan and Mayr, 1989]. This type of 
ionopause layer was observed occasionally during solar maximum 
and is typical during solar minimum. 

While it is generally agreed that the ionopause is a layer of finite 
thickness across which the plasma changes from that originating 
predominantly in the sun to that originating predominantly in the 
planet Venus, some experimenters have been more specific in 
defining the ionopause as a surface which lies within theionopause 
layer. For example, Knudsen et al. [ 1982b 1 used the term, top of 
the ionosphere, to refer to the bottom of the ionopause layer for low 
solar wind dynamic pressure transitions and the term, top of the 
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Fig. 4. Altitude profiles of the major ions 0 + and 0 2 * measured on the 
outbound legs of two orbits. The normal component of the solar wind 
dynamic pressure, obtained from magnetometer measurements, was low, 
1.8 nPa, for orbit 420 and high, 8.5 nPa, for orbit 412 (adapted from 
[Mahajan and Mayr, 1989]). 

ledge to refer to the top of the layer. Brace et al. [1980] have 
defined the ionopause as the altitude where the plasma density 
decreases below approximately 100 cm' 3 . This location is 
approximately equivalent to Knudsen’s top of the ledge. Elphic et 
al. [1980] and Phillips et al. [1988] define the ionopause altitude as 
that at which the magnetic pressure in the overlying magnetosheath 
first equals the thermal pressure of the ionospheric plasma. 

The dynamic pressure of the solar wind is transmitted to the 
ionosphere primarily through the pressure of the magnetic barrier. 
As the shocked solar wind approaches the ionopause, the embedded 
magnetic field lines become draped around the ionopause and lose 
much of their plasma [Zwan and Wolf, 1976], The region of large 
magnetic pressure and small particle pressure at the top of the 
ionopause (top of the ledge) is called the magnetic barrier. Figure 5 
illustrates the variation of magnetic and ionospheric thermal 
pressure during one periapsis pass of the PVO spacecraft through 
the ionosphere. At an altitude above 400 km, the pressure is 
essentially all magnetic. At the bottom of the transition layer, the 
magnetic pressure has become negligible when compared to the 
ionospheric thermal pressure. At the lowest altitudes in this 
example, the magnetic field pressure again builds up to about 30% 
of the total, hut this field plays no important role in the .structure of 
the ionopause layer. Figure 5 illustrates the behavior of the 
magnetic field and ionospheric plasma in the ionopause transition 
layer for cases of moderate to low solar wind dynamic pressure. 
For cases with high solar wind dynamic pressure, the magnetic field 
remains large within the ionosphere down to the altitude of 
periapsis and can affect the ionospheric plasma [Luhmann et al., 
1980]. 

The variation of the "pressure balance" ionopause altitude with 
ionopause pressure is illustrated in Figure 6. Total ionopause 



pressure for this figure is the sum of magnetic pressure and 
ionospheric thermal pressure measured at the ionopause height. 
The total pressure is presumed to be equal to the component of the 
solar wind dynamic pressure normal to the ionopause surface 
reduced by the factor K (= 0.86). For total pressure less that about 
4 nPa, the median ionopause altitude decreases approximately 
linearly with increase in the logarithm of the total pressure. This is 
the behavior expected if the ionospheric particle pressure is 
balancing the total ionopause pressure. The logarithm of the 
ionospheric thermal pressure decreases approximately linearly with 
altitude (Figure 7). As the total pressure increases above 4nPa, the 
median ionopause altitude approaches a lower limit of 
approximately 250 km. The large solar wind dynamic pressure 
case illustrated in Figure 3 is an example of an ionopause at the low 
altitude limit. For total ionopause pressures of approximately 4 nPa 
and larger, the ionospheric particle pressure is evidently unable by 


itself to balance the total ionopause pressure. 
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Fig. 6. Variation or ihc ionopause height wiih total ionopause pressure at 
solar maximum (adapted from (Phillips cl al., 1988J). Total ionopause 
pressure is the sum of magnetic and ionospheric thermal pressure at the 
height of equality. All crossings with SZA less than 120° arc plotted. Solid 
diamonds are median altitudes of log pressure bins. 
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Fig. 7. Median profiles of the ionospheric particle pressure n k(T +T > 
during solar maximum. Medians are medians of ,he local product and no 
he product of .he medians. Note .ha. .he altitude scale is different for .he 
two solar zenith angle internals (from [Miller et a!., 1984]). 


The large variation in ionopause heights for a given total 
pressure illustrated in Figure 6 suggests that dynamical processes, 
such as the Kelvin-Helmholtz instability and other as yet 
unrecognized processes arc operating lo produce instantaneous 
ionopause heights that differ from that expected under steady stale 
conditions. 

Figure 8 illustrates the relative frequency or occurrence of the 
up-stream solar wind dynamic pressure during 1980, a period of 
solar maximum and the "effective" peak ionospheric pressure 
(6x10 dyne/cm ) for the same period. The "effective" peak 
ionospheric pressure is the median peak ionospheric panicle 
pressure of Figure 7 divided by the factor K. Figure 8 illustrates 
that the up-stream solar wind dynamic pressure was less than the 
effective peak ionospheric particle pressure about 80% of the time 
during solar maximum. 


- - r — — •vMvzpaujt uuiing boiar maximum can 
be accurately modeled by assuming that the ionopause altitude is 
determined by pressure balance between the sum of the solar wind 
dynamic and thermal pressures and the ionospheric thermal 
pressure. Figure 9 illustrates the computed altitude of the 
ionopause location and the measured median altitude at solar 
maximum for SZA between 0°and 125° [Knudsen et al 198^1 
The ionospheric pressure Held was derived from measured 
quantises. The computed and measured median ionopause heights 
agree well at solar maximum because the median up-stream solar 
wind dynamic pressure is moderate, and the ionospheric thermal 
pressures able to stand off the solar wind most of the time 
(Figure 7). 


Phillips et al. f 1988] have recently shown that a small cylindrical 
no™*'- 0 ' re P° rted b y Br ^e et al. [1980] and Knudsen et al. 
[1982b] in which the dawn ionopause altitude is higher than the 
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dusk ionopause is largely the result of orbital bias and uneven 
sampling of solar wind dynamic pressure. However, they found an 
asymmetry in which the ionopause in the terminator region is 
displaced in a direction opposite to that of the motional electric 
field. 

The solar maximum ionopause shape beyond about 125° SZA 
has not been determined experimentally because of measurement 
limitations imposed by the PVO orbit. However, using a gas 
dynamic numerical model, Moore et al. [1991] have recently 
derived an ionopause shape extending to a large distance down the 
Venus wake. Figure 10 illustrates the shape of the ionopause that 
he obtained by iteratively altering the.shape until the solar wind 
normal pressure at the ionopause balanced the ionospheric pressure 
for x between 0 and -0.5. 

Wave-like structure observed within but near the ionopause, as 
illustrated in the electron density and temperature and magnetic 
field for orbit 184 of Figure 1J, has been interpreted as evidence for 
waves on the ionopause [Brace et al., 1980, 1982a]. In some 
instances a "cloud” of plasma, such as that occurring on the inbound 
leg of orbit 173, appears to be completely detached from the main 
body of the ionopause. The strong decrease in density, increase in 
magnetic field strength, and increase in T c all suggest that the 
spacecraft has reentered the ionosheath after encountering a patch 
of ionospheric plasma. Some patches of plasma are encountered 
with a greater separation between the cloud and the ionopause than 
is indicated in orbit 173. Although there is separation of the cloud 
from the ionopause along the orbital tract, the clouds may not be 
completely detached from the ionosphere proper. The orbital track 


is near-tangent to the ionopause and irregularities in ionopause 
height could also appear as detached plasma. 

The OIMS has observed supertherma! ions near and above the 
ionopause as illustrated in Figure 12. The behavior of the 
instrument current in the shaded region above 270 km altitude is 
interpreted as consisting of ions with a kinetic energy component 
parallel to the instrument axis in the range 30-80 eV. The inference 
is that 0 + ions in the shaded region have been accelerated 
antisunward into this energy range through interaction with the 
rapidly flowing solar wind plasma [Taylor et al., 1980a]. 

Spenner et al. [1980] found that the differential electron energy 
flux spectrum during solar maximum changed from one 
characteristic of the upper ionosphere photoelectron spectrum to 
one characteristic of the ionosheath solar wind electrons as the PVO 
spacecraft passed upward through the ionopause and into the 
ionosheath. Figure 13 illustrates spectra obtained in the ionosheath, 
in the ionosphere and in the transition layer. They called this layer 
the mantle (see Figure 1) and suggested that the upper surface of 
the mantle was the effective boundary diverting the shocked solar 
wind flow around the planet. 

Density and Composition 

Figure 14 illustrates the altitude and solar zenith angle 
dependence of the median ionospheric total ion density during solar 
maximum. The density varies approximately as the square root of 
SZA over most of the dayside hemisphere but drops rapidly by 
about an order of magnitude across the the terminator. Figure 15 
illustrates the order of magnitude drop in electron density (equal to 
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Fig. 10. Ionopause shape extrapolated into the near wake. The 
extrapolation is based on pressure balance between the ionospheric thermal 
pressure and solar wind thermal pressure over the region 0 > x > -0.5 where 
x is measured in units of the Venus radius [Moore et al., 1991). 

total ion density) across the terminator at three discrete altitudes. 
The dots are individual density measurements made during the first 
600 orbits within the altitude intervals 160-170, 260-280, and 550- 
650 km. They reveal the strong increase in temporal density 
variability from the dayside to nightside hemispheres. The solid 
lines are empirical model values of the density [Theis et al., 1984], 

The solar maximum dayside ionosphere of Venus consists 
principally of an upper layer of 0 + ions and a lower layer of 0 2 + 
ions. Other ions compose less than 10% of the total. Figure 16 
illustrates the major ion composition measured by the ORPA for an 
SZA of 65°. The vertical profiles are smoothed median density 
profiles and the temporal variation bars indicate the interval within 
which two thirds of the measured values fall. Only measurements 
lying below the ionopause were included in the data base. The ion 
density labeled M29 + density is the sum of the densities of ions with 
mass 28 amu (N 2 + + CO + ) and 30 amu (NO + ). The ORPA does not 
resolve M29 + from 0 2 + , and the density of these two individual 
ions is less certain than their sum. The 0 + , 0 2 *, and C0 2 + densities 
vary as cos(SZA) 1 ^ 2 across most of the dayside hemisphere [Miller 
et al,, 1984]. Figure 17 illustrates the ion composition measured by 
the PVO OIMS during one subsolar orbit passage through periapsts 
[Taylor et al., 1980a]. The ion labeled H 2 + in Figure 17a is now 
known tobeD + [McElroy et al., 1982; Donahue et al., 1982; Hartle 
and Taylor, 1983]. The solar zenith angle variation of several jon 
species at an altitude of 200 km measured by the OIMS is presented 
in Figure 18. Perhaps the most interesting variation is the strong 
increase in H + ion density. The H + ion density becomes comparable 
to 0 + in a range of SZA centered on 130° in the dawn hemisphere 
in consequence of the enhanced neutral hydrogen density at this 
location and charge exchange between 0 + and H [Brinton et al., 
1980]. 

During solar maximum, the nightside ionosphere, like the 
dayside ionosphere, typically consists of a layer of 0 + ions 
overlying a lower Of layer as illustrated by the median vertical 
profiles in Figures 19 and 20. The variance bars reflect the fact that 
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Fig. 11. (a) Measurements of N c , T c , and B from orbit 173 showing the 
passage of the PVO spacecraft back into the jonnshealh (magnetic harrier ?) 
for a short period of time after first encountering the ionosphere as the 
spacecraft entered the ionopause at a shallow angle, (b) Measurements of 
N c , T c , and B exhibiting wave like character as the spacecraft transited the 
ionopause inbound and outbound on orbit 184 (from [Brace et al.. 1980]). 

the temporal variability of the nightside ionospheric density is 
larger than that of the dayside ionospheric density (compare with 
Figure 15 and 16). In the SZA range 120°-150° of the dawn 
hemisphere, the H + density becomes comparable to and greater than 
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Fig. 12. A daysidc example (SZA = 11 °) of the plasma measurements from 
the OIMS, including the bowshock-ionosheath, with energetic ion currents 
(energy > lOOeV) peaking near 1700 km, a region of superthermal 
(energy = 30-80 eV) ion flow (shaded) extending across the range 270-430 
km, and the relatively undisturbed thermal ionosphere in the interval from 
270 km down to periapsis at 173 km. The ionopause is identified at 
approximately 270 km, the altitude above which the superthermal ions are 
observed. Within the shaded region thermal ion concentrations cannot be 
precisely determined owing to incomplete sampling of three dimensional 
plasma flow characteristics. The shaded area indicates qualitatively the 
height dimension and relative intensity (not calibrated in either density of 
species) of the superthermal flow (from [Taylor et a!., 1980a]). 



Fig. 14. Global variation of the median ionospheric ion density at solar 
maximum (after [Miller et al., 1980]). 


the 0 + density as shown in Figure 18a for the reasons presented in 
the previous paragraph. A vertical profile of the median Ff density 
in the region of enhancement is illustrated in Figure 21. In 
comparing the median H* vertical profile of Figure 21 with that of 
0 + in Figure 19, note the difference in altitude scales of the two 
figures. 

Figure 22 illustrates the ion composition measured by the OIMS 
along orbit 65 outbound from its periapsis at SZA = 162°. The 
composition on the nightside is, for the most part, similar to that on 
the dayside, but the density level is about a factor of ten less. The 
significantly higher density of FT is an exception. The similarity in 
composition, and the observation of large ion bulk velocity directed 
antisunward at the terminator led several of the PVO experimenters 



Fig. 13. Change in differential electron flux as the PVO spacecraft passed from the ionosphere into the ionosheath (magnelosheath) 
at solar maximum (from [Spenner et al., 19801). The spectrum changed from one characteristic of ionospheric photoelectrons to one 
characteristic of the ionosheath electrons. 
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Fig. 15. Total electron density measured by the OETP within three altitude 
intervals, 160-170, 260 - 280, and 550 - 650 km as a function of SZA. The 
solid lines are empirical model values at the indicated altitudes derived from 
median values of the individual measurements, the latter indicated by the 
dots. The data base consists of measurements within the ionosphere from 
Dec 1978 to Dec 1982, the period of solar maximum (from [Theis et al., 
1984]). 



ION DENSITY (cm" 1 ] 

Fig. 16. Vertical profiles of the median density of the major ions measured 
in the dayside hemisphere by the PVO ORPA during solar maximum. M29 + 
represents the sum density of ions with masses 28 and 30 amu. Together, 
0 + and 0 2 + compose in excess of 90% of the density (adapted from [Miller 
et al., 1984]). 


fa 1 


PIONEER VENUS OIMS ORBIT 185 SZA-IT DAY 



10° lO 1 IO* 10 s I0 4 I0 5 I0 6 

ION DENSITY (tONS/CM 5 ) 


lb) 

PIONEER VENUS OIMS ORBIT 105 SZA-II - DAY 



K>' 10* I0 3 I0 4 I0 5 I0 6 


tON DENSITY (I0NS/CM 3 ) 

Fig. 17. A dayside example of the major and minor ions measured in the 
lower ionosphere of Venus by the PVO OIMS at solar maximum ([Taylor 
et al., 1980a]). Outbound pass occurred near local noon with SZA = 11° at 
the location of periapsis. Distributions are divided into panels (a) and (b) 
for clarity. 

to postulate that transport of dayside ionospheric plasma into the 
nightside ionosphere was a significant source of the nightside 
ionosphere (Taylor et al., 1979a; Knudsen et al., 1980b], 

The density and composition of the dayside ionosphere sunward 
of the terminator is reproduced reasonably accurately in one 
dimensional numerical models by assuming ion production from 
solar photon and photoclectron ionization of the neutral atmosphere 
and loss through a chain of chemical reactions. In modeling the 
behavior of ions above 170 km altitude, it is necessary to include 
vertical diffusion fKnudsen et al., 1979a; Nagy et al., 1979; Nagy 
et al., 1980]. At solar maximum, horizontal transport of ionization 
across the terminator is essentia! in modeling the nightside 
ionosphere, especially the 0 + layer, and several two dimensional 
numerical models have been developed to predict the major ion 
densities as a function of SZA [Whitten et al., 1982; Cravens et al., 
1983; Whitten et al., 1984; McCormich et al., 1987; Slnghal and 
Whitten, 1987]. These authors recognized that electron impact 
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contributes to the nightsidc ionization, but this source was not 
included in the models. 

Energetics 

The global variation of electron temperature (T e ) measured by 
the OETP at selected altitudes is presented in Figure 23. The dots 
are all the individual measurements from the first 600 orbits, a 
period of solar maximum, and the solid lines are from an empirical 
model derived from the medians of a subset of the data base. For 
comparison. Figure 24 illustrates a vertical profile of median values 


of T c measured by the ORPA in several SZA ranges. There is 
evidently little variation of the median T e with SZA over the entire 
planet, although Figure 23 does exhibit a small decrease from the 
dayside to the nightside hemisphere. Numerical model calculations 
of T e over the dayside hemisphere have shown that the measured T e 
profile cannot be reproduced when solar EUV heating is the only 
heat source unless the thermal heat conductivity is greatly reduced 
by the presence of steady and fluctuating magnetic fields [Cravens 
et al. f 1980, Kim et al„ 1990]. By assuming a heat inflow at the 
ionopause of the order of a few times 10 10 eV cm'\s' ( it is possible 
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Fig. 19, Vertical profile of median 0 + density in the solar maximum 
nightside ionosphere and mode! computations of the 0 + layer expected from 
a downward flux of 0 + ions and from electron precipitation (from [Knudsen 
et al„ 1986]). 


to match the observed temperature profile [Knudsen et al., 1979b; 
Cravens et al„ 1980; Kim et al., 1990], 

Figure 23 illustrates the large temporal variability of T e , 
especially in the night hemisphere. A significant part of this 
variability results from inclusions in the data base of T e 
measurements when the electron density was very low, such as in 
holes and during periods of greatly depressed ionospheric density 
(disappearing ionosphere), topics to be discussed later. 

The ion temperature Tj shows little variation with SZA over the 
dayside hemisphere, but changes substantially near the terminator 
and in the nightside hemisphere (Figure 25). The variation between 
170 and 200 km altitude on the dayside requires the existence of a 
heat source localized in altitude and restricted to the dayside 
hemisphere. It has been identified as that due to chemical heating 
[Rohrbaugh et al., 1979; Cravens el al., 1979; Knudsen et a!., 
1979b; Kim et al., 1990]. Reproduction of observed dayside ion 
temperatures above 200 km altitude with numerical models 
requires a distributed heat source or a heat flux at the ionopause 
unless the ion heat conductivity is reduced by a horizontal and/or 
fluctuating magnetic field [Knudsen et al., 1979b; Cravens et al., 
1980;, Kim et al., 1990], The increase of Tj with SZA in the altitude 
range 170-300 km at and beyond the terminator is probably due to 
decreased ion heat loss to the neutrals. (The neutral density and 
temperature decrease rapidly with SZA across the terminator.) The 
high altitude ionosphere flows supersonically across the terminator, 
as we shall discuss later,^md advection of heat within the ion gas 
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Fig. 20. Vertical profile of the median 0 2 + + M29* sum density in the solar 
maximum nightside ionosphere. The symbol M29 + is defined in the 
Figure 16 caption (from [Knudsen et al., 1986]). 
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Fig. 21. Vertical profile of median H + density in solar maximum nightside 
dawn hemisphere between 120 and 150 SZA (from [Knudsen et al., 1986)). 
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Fig. 22. A nightside (dusk hemisphere) example of the ion composition 
measured by the OIMS in the outbound leg of orbit 65. The SZA at 
periapsis was 162°. The strong decrease in density at 230 km altitude may 
result from passage through an ionospheric hole (from [Taylor et al., 
1980a]). 



Fig. 23. Total electron temperature measured by the OETP within three 
altitude intervals, 160 -170, 260 - 280, and 550 - 650 km as a function ol 
SZA. The solid lines are empirical model values at the indicated altitudes 
derived from median values of the individual measurements, the latter 
indicated by the dots. The data base consists of measurements within the 
ionosphere from Dec 1978 to Dec 1982, the period of solar maximum 
(adapted from [Theis ct al., 1984]). 



ELECTRON TEMPERATURE (K) 


Fig. 24. Vertical profiles of median electron temperatures as a function of 
SZA at solar maximum (from [Miller et al., 1980]). 


and conversion of the bulk kinetic energy into heat become 
important [Knudsen et al M 1981a; Bougher and Cravens, 1984; 
Whitten et a!., 1986]. 

Velocity Field 

Rapid advection of the near-terminator dayside 0 + layer into the 
nightside hemisphere plays a major role in the formation of the 
nightside ionosphere and in the nightside ion energetics. Figure 26 
illustrates the median velocity field of 0 + ions measured on the 
outbound portion of orbits from periapsis to the ionopause. These 
data represent near-equatorial median results obtained from all 
solar hour angles over a period of approximately three and a half 
Venus years [Knudsen el al., 1982a]. The bulk velocity of the ions 
becomes reasonably well organized near the terminator and has a 
magnitude of the order of 3-4 km/s at higher altitude and 
antisunward of the terminator. The median velocity vectors are 
generally smaller and confused in the sub- and antisolar regions in 
this presentation. 

Some of the confusion in the subsolar region disappears when 
the data are separated into dawn and dusk hemispheres before 
taking averages as shown in Figure 27. The altitude above the 
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Fig. 25. Vertical profiles of median ion temperatures as a function of SZA 
at solar maximum (from [Miller et al, 1980]). 


planet is exaggerated by a factor of four relative to the radius of the 
planet. In this analysis the flow in the subsolar region and at lower 
altitudes is seen to be westward and in the direction of planet 
rotation. The rotation of the cloud tops is in this direction and of 
the order of 0.1 km/s. If the thermosphere were super rotating at a 
speed of the order 1 km/s, we would expect this observed behavior. 
Because no corroborating evidence or theory for such a large super 
rotation exists, this result must be considered tentative, * An 
uncertainty sphere of 0.3 km/s radius exists at the head of 
individually measured ORPA velocity vectors under the best of 
circumstances and some as yet undetected systematic error may 
also be present. 

The measured average velocity vectors become chaotic beyond 
about 150° SZA in both Figures 26 and 27, The reason for this is 
not established at the moment but several factors may be 
contributing. The number of measurements at the higher altitudes 
with SZA greater than 150° is smaller than at lower SZA because of 
orbital coverage and also because of failure of the ORPA to detect 
O 4 ion peaks at the low densities and high ion temperatures 
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Fig. 26. Bulk ionospheric velocity field at solar maximum (from (Knudscn 
et al., 1982a]). 


occurring in that region. The horizontal component of any 
cylindrical^ symmetric flow field about the sun-Venus axis should 
approach zero in the antisolar region. It is likely that the 
transterminator flow field is not cylindrically symmetric near the 
terminator and that the center of convergence of the flow wanders 
around in time. All these factors are probably contributing to the 
confused average flow near the antisolar axis. 

Theis et al. [1984], Whitten et al. [1984] and Elphic et al. [1984], 
using two dimensional numerical models, have been able to 
reproduce reasonably well the measured velocity field. All models 
indicate that the measured pressure field is able to accelerate the 
plasma to the observed velocities. In the model by Whitten et al. 
[1984] the pressure field is produced entirely by EUV ionization. 
This result indicates that momentum transfer from the solar wind 
via viscous drag is not important in accelerating the ionospheric 
plasma [Knudsen et al., 1981b]. 

For flow speeds in excess of about I km/s, the flow is supersonic 
| Knudsen et al., 1980b]. Knudscn ct al. [ 1980b] have suggested lhal 
a shock wave forms in the vicinity of 150° SZA which slows ihc 
flow, changes its direction and heats the ions in the central nightside 
region (Figure 28). This suggestion has recently been supported by 
observations of median ion temperature and density increases of the 
expected magnitude and at the expected SZA, and by the results of 
a model [ Miller el al., 1990; Kim et al., 1990J. 

Nightside Ionospheric Holes 

The term ionospheric hole has been given to an intriguing and 
unanticipated phenomenon discovered in the nightside Venus 
ionosphere by Brace et al. [1982b]. They observed that strong 
plasma density depletions frequently occurred within small time (or 
space) intervals during nightside periapsis passes. Outside these 
time intervals the plasma density was at its undisturbed value. 
Figure 29 illustrates the plasma density measured during several 
passes through the nightside ionosphere. The ionospheric holes are 
evident in several of these orbits. Two holes, one on either side of 
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Fig. 27. Bulk ionospheric velocity field separated into dawn and dusk hemispheres. The altitude scale above the planet surface is 
four times that of the planet. The low altitude field in the dawn hemisphere suggests that a neutral atmospheric super rotation may 
be exerting an influence on the flow (from [Miller and Knudsen, 1987]). 
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Fig. 29. Ionospheric holes present in several orbits that otherwise exhibit a 
fully established nightside ionosphere (from [Brace et al„ [982b]). 

periapsis, ;irv evident in each of the orbits 530, 531, 533 and 537. 
The holes are not seen in every passage and appear to move in 
latitude. The density depression tends to be greater in the holes 
occurring in the outbound portion of pcriapsis passes than those in 
the inbound portion which suggests that the holes disappear at 
altitudes near and below the nightside main peak altitude. 
Figure 30 illustrates the presence of an enhanced magnetic field 
within the holes directed predominantly along the sun-Venus axis. 

Marubashi et al. [1985] have projected the location and polarity 
of the magnetic field observed in 55 holes onto the nightside planet 


after rotation into a solar wind magnetic coordinate system so that 
the upstream Y component of the solar wind magnetic field is 
horizontal (Figure 31). Holes with sunward directed magnetic 
fields are plotted with crosses and holes with antisunward pointing 
fields with circles. This rotation orders the holes with sun-directed 
magnetic fields on the right and those with antisunward directed 
fields to the left of the Y=0 plane. This ordering is consistent with 
the fields having their origin in the draped magnetosheath field as 
indicated on the right. Processes responsible for the formation of 
the holes have been suggested by Grebowsky and Curtis [1981], 
Brace et al. [1982a], Luhmann et al. [1982], Grebowsky et al. 
[1983] and Marubashi et al. [1985]. 

Post Terminator Waves 

Brace et al. [1983b] have reported the presence of wave-like 
structures in measurements of electron density, electron 
temperature, and magnetic field at low altitudes antisunward of the 
terminator. Figure 32 illustrates the waves and the anti-correlation 
between T c and N e , and correlation of peaks in N c with sign 
changes in B E , where the magnetic field components are radial 
eastward B E and northward B N . These waves first appear at the 
terminator, extend for at least another 30° SZA into the night 
hemisphere, and exist primarily below 200 km altitude. The waves 
are evidently not accompanied by waves in the neutral atmosphere 
[Brace et al., 1983b]. The generation mechanism for the waves is 
unknown, but the fact that they occur in a region of large plasma 
velocity shear suggest that a shear instability may be involved. 

Suprathermal Electrons 

Although the PVO spacecraft was not adequately instrumented 
to measure energetic particles within the ionosphere with energies 
above about 1 eV energy, the low energy ionospheric and neutral 
atmospheric instruments have been able to provide some 
measurements of considerable value for understanding the Venus 
energetic particle environment. 

The ORPA has a mode of operation in which electron current is 
collected at several retarding voltages from 0 to 48 V. The 
differential photoelectron flux derived for the dayside solar 
maximum ionosphere from such measurements is illustrated in 
Figure 33. Comparison of the differential photoelectron fluxes 
computed numerically [Cravens et al., 1980] with the measured 
flux suggests that the measured flux is produced predominantly by 
solar photons and that the photoelectrons diffuse rather freely in the 
vertical direction. 

Spenner et al. [1981] and Knudsen and Miller [1985] have 
reported on the suprathermal electron (lux observed by the ORPA 
in the Venus nightside ionosphere during ihe l f )K() 19X1 snlar 
maximum and the following waning phase of the solar cycle. The 
measurements used in the study were restricted to the altitude 
interval of l(XX)-20(X) km and to SZA greater than I3(P. a region 
similar to that sampled by the Venera 9 and 10 suprathermal 
electron detectors during an solar minimum period [Gringauz cl al., 
1979]. The Venera instrument had a retarding potential range of 
0-300 V. Figure 34 illustrates the median electron spectrum derived 
in the study by Knudsen and Miller and compares it with spectra 
reported by Gringauz et al. [1979] for the same region but measured 
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Fig. 30. Enhanced radial magnetic field in ionospheric holes. The magnetic field components are in solar ecliptic coordinate system; 
the angle V is between the radial direction and the local magnetic vector. M 7 = 0 is radial. The magnetic field is enhanced within the 
holes and is oriented at a small angle to the radial. The vertical lines delineate the boundaries of the holes (from (Luhmann et al., 

1982]). 

during the preceding solar minimum period. Both Knudsen and depended on the assumed temperature (T 3 ). For an assumed T 3 
Miller and Gringauz et al. found that the supratherma! flux could be temperature of 200 eV the density was 0.06 cm 3 . The median 

approximated by a two component Maxwellian electron gas. values of n 2 , T 2 and T 3 were found to be uniform over the altitude 

Knudsen and Miller found that the most abundant component has a and SZA ranges covered by the study. Two thirds of the measured 

median density (n 2 ) of 3 cm " 3 and a median temperature (T 2 ) of values of n 2 , T 2 and n 3 fell within a factor of +/-2 of the median 

14 eV. The density of the high temperature component (r^) values. Comparing the ORPA median omnidirectional flux with 
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Fig. 31. PVO orbit segments where holes were detected are projected onto the nightside planet after rotation into the Y-Z plane of 
the solar wind magnetic coordinate system (left). Holes with magnetic field directed sunward are plotted with crosses, and holes 
with magnetic field directed antisunward, with circles. This organization of the segments is consistent with the hole magnetic fields 
originating in the draped solar wind magnetic field (right) (adapted from [Mambashi et al. t 1985]). 
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Fig. 32. The ionospheric and magnetic characteristics of post-terminator 
waves in orbit 473 (SZA=100° at periapsis). The upper five traces show a 
6-minute segment of N c , T c , B R , Bp, B y , and iBi, while the lower panels 
display and expanded view of N c . T c , and Bp which show details of the wave 
patterns. Well ordered waves develop below about 175 km altitude. The 
vortical dashed lines draw attention to persistent phase relationships 
between Bp. N c . and T c (from |Brace el :il„ l«)83h|). 


some statistical data reported by Gringauz et a!., Knudscn and 
Miller concluded that the ORPA median flux was about a factor of 
2 larger than the Venera-measured fluxes within the lower energy 
range common to the two instruments, in both studies the 
omnidirectional flux was derived by assuming the flux to be 
isotropic. 

Little directional information on the suprathermal flux is 
available. The effective solid angle of the ORPA is about pi 


steradian, and its axis lies on a cone with a half angle of 25°. The 
axis of the cone is perpendicular to the earth ecliptic. The density 
and temperature of the suprathermal electrons are derived with the 
assumption that the electron flux is isotropic. No statistical studies 
have been conducted at the present of the ORPA suprathermal 
electron flux to look for a possible solar cycle variation or for flux 
anisotropies. The Venera electron detector has about the same solid 
angle as that of the ORPA, but the axis was fixed in the antisunward 
direction. Thus, all that can be inferred at the present is that the flux 
is probably approximately isotropic and that there is not a large 
solar cycle variation. 

Nightside Ionization Sources 

Two ionization sources play a significant role in maintaining the 
nightside ionosphere. The first to be recognized was the electron 
impact source, a source produced by the aforementioned 
suprathermal electron fluxes with energy of the order of 14 eV. 
Gringauz et al. [1979] were the first to recognize that the 
suprathermal electrons were a source of ionization and that the 
measured flux could produce the nightside electron density profiles 
measured by the Venera radio occultation experiments provided the 
electrons impacted a neutral CO 2 layer with appropriate scale 
height and density. The suprathermal electron fluxes are present in 
the nightside Venus hemisphere at solar maximum and solar 
minimum and contribute, therefore, to the maintenance of the 
nightside ionosphere at all phases of the solar cycle. The PVO 
measurements suggest that the median electron impact source may 



Fig. 33. Differential photoclcclron flux in ihe solar maximum dayside 
ionosphere (from (Knudscn et a!., 1980aj). 
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Fig. 34. Median supralhermal electron energy spectra measured in the 
Venus umbra at solar maximum in the altitude interval 1000-2000 km by the 
PVO ORPA and at solar minimum by the Venera 9 and 10 hemispherical 
retarding potential analyzer ( adapted from [Knudsen and Miller, 19851). 

be larger by a factor of two at solar maximum than at solar 
minimum. Gringauz et al. and Knudsen and Miller [1985] have 
both reported considerable variability in the suprathermal electron 
fluxes. The variability, observed within single orbits as well as 
from orbit to orbit, may be reasonably assumed to result from 
temporal as well as spatial variations. Recent model computations 
of the ionization produced by the measured suprathermal electron 
fluxes impacting the solar maximum neutral nightside atmosphere 
demonstrate that these fluxes are primarily effective in producing 
an 02 + ionospheric layer [Spenner et al., 1981; Cravens et al., 
1983]. The O* layer peak electron density produced by the flux is 
smaller by a factor of 10 2 than the 0 2 * la y er Pe ak electron density 
(see Figures 19 and 20). The computed column ion production rate 
from electron impact in both studies was ol the order ol 
UI0 7 u»’js 

The second ionization source to be recognized was the () 
transport source. Taylor et al. [ 1979|, observing that the nightside 
composition was similar to that in the dayside and that large 
horizontal ion velocities existed, concluded that transport of CT ions 
could produce the observed nightside ionization. The Of main 
peak would be produced through the reaction 
0 + + C0 2 -*0 2 + +C0. Knudsen et al. [1980b] established the 
importance of 0 + transport by making an ionization balance 
calculation. From early ORPA vector velocity and 0 + density 


measurements, they calculated the median flux of 0 + ions crossing 
the terminator to be 8x10 26 s* ] (Table I). The median column 
source flux Sj that this transterminator flux would provide if 
distributed evenly across the the nightside hemisphere is 
3 xl0 8 cm' 2 s’ 1 . They also estimated the average column 
recombination rate L across the nightside hemisphere from the 
shape of the 0 2 * layer s(r) and its peak density n^G) as a function 
of SZA to be 2xI0 8 cm* 2 s' 1 . Thus, the transterminator 0 + flux w as 
sufficient to maintain the nightside hemisphere provided that most 
of the flux was being distributed appropriately across the entire 
nightside hemisphere and returning to the planet. Later studies have 
shown that the 0 + flow typically canies to at least 135° SZA (see 
Figure 26) and, hence, covers the bulk of the nightside hemisphere. 
Present evidence suggests that the planetary ionization escaping 
down the Venus wake is smaller by more than an order of 
magnitude than the cross terminator 0 + flux [Knudsen et al., 1980b, 
Mihalov and Barnes, 1982; and McComas et al., 1986]. 
Consequently, most of the transterminator 0 + flux recombines in 
the ionosphere. Numerical model studies, one dimensional as well 
as two dimensional, have demonstrated that formation of the 
nightside 0 + layer requires a downward flux of 0 + ions with 
magnitude of the order of IxlO 8 cm' 2 s’*, and that the flow w ill 
carry through most of the nightside hemisphere [Spenner et al., 
1981; Cravens et al., 1983; Whitten et al., 1984; McCormick et al.. 
1986; and Nagy et al., 1991]. 

With the discovery of the 0 + transport source, the question arose 
as to the relative significance of the two sources. Spenner et al. 
[1981] concluded that the 0 + transport source was dominant with 
the 0 + transport source typically supplying three times the 
ionization that electron impact was supplying. The evidence cited 
above suggests that this conclusion was conservative for typical 
conditions existing during the first few years of the PVO mission. 
However, evidence is presented hereafter which shows that the O 
transport source is typically large only during solar maximum and 
is probably insignificant during solar minimum. The two sources 
evidently alternate in being the dominant nightside ionization 
source. 


TABLE I. Nightside lonoshperic Ionization Balance 


F , 
(I0 2fi s* 1 

Si s e 

) (loWV) (lO^cmV) 

L 

(10 R cm' 2 s 

Reference 

;•') 

8 

3 

2 

Knudsen et al. 1980 


-.02 

2 

Spenner et al. 1980 




Cravens et al 1983 


;ii«mI lumsicnninnior <>’ Hux “ n 2 fti Ar 


S ( : Average column source = 17(2nK/) 

S c : Column electron impact ionization rale 
L: Nightside average column recombination rate 

L = —L [" [~ an 2 (6, r) 2nr 2 sin0d6dr 

2ltR 2jn/2jRv 


n( 0 , r) =n m (0)s(r) 






256 The Venus Ionosphere 


Night side Aurora 

The Pioneer Venus ultraviolet spectrometer has observed 130.4 nm 
radiation, spatially and temporally variable in brightness and in 
coverage, emanating from the nightside hemisphere [Phillips et aJ. t 
1986], The variability and energy spectrum of the suprathermal 
electrons observed in the nightside ionosphere are consistent with 
the suprathermal electrons being the source of the aurora [Fox and 
Bougher, 1991; Fox and Stewart, 1991]. No direct correlation of 
the 130.4 nm brightness with suprathermal electron flux is possible 
since the EUV observations are made when periapsis is in the 
dayside hemisphere and the suprathermal electron flux 
measurements when periapsis is in the nightside hemisphere. 

Nightside Density Variability 

After the first few passages of the PVO spacecraft periapsis 
through the night hemisphere, a period of solar maximum, it 
became evident that the ionosphere was substantial or "established” 
most of the time. However, it was almost unobservable at the PVO 
altitude on some occasions and intermediate on other occasions 
[Taylor et al., 1979; Knudsen et al, 1979b], Figure 29, introduced 
earlier, illustrates the variation of N e with time as the PVO 
spacecraft passed through the nightside hemisphere during periods 
when the ionosphere was established. For comparison. Figure 35 
shows an extreme example of what has been called a 
"disappearing" nightside ionosphere. The greater variability of the 
nightside ionospheric density relative to that in the dayside 
hemisphere is evident in the scatter of points in Figure 15 and in the 
sizes of the temporal variance bars in Figures 16, 19 and 20. 

The cause of the apparent nightside variability has to be 
variability in the two sources of ionization or in the loss processes. 
Gringauz et al. [1979] concluded that the variability in the nightside 
Venera radio occultation profiles, which were recorded during solar 
minimum, was caused by the observed variability of suprathermal 
electron fluxes. Since we now have evidence that the CF transport 
ionization source is greatly reduced at solar minimum, their 
conclusion seems reasonable. 

Knudsen et al. [1980b] and Spenner et al. [1981], recognizing 
that transterminator flow of ions from the dayside ionosphere was 
a dominant source of the nightside ionosphere at solar maximum, 
suggested that the relatively large variability of the nightside 
ionosphere was caused by temporal and spatial variations in the 
transterminator 0 + flux. Since that flux depends on the ionopause 
height al the terminator, as well as on the velocity field, it is 
reasonable to assume that a correlation might exist between the 
density of the nightside ionosphere and the terminator ionopause 
height. Two dimensional model studies in which the height of the 
ionopause at the terminator has been varied confirmed this 
expectation [Cravens et al., 1983; Knudsen et al., 1987], Miller and 
Knudsen [1987] investigated the correlation between ionopause 
height and nightside electron density experimentally. In Figure 36, 
the ionopause height measured near the terminator and normalized 
by the median ionopause at the same SZA is plotted against the ion 
density measured in the same orbit at an altitude between 170 and 
200 km altitude and normalized by the local electron density of the 
empirical model of Theis et al. [1984], The selected orbits had 
penapses in the central nightside hemisphere. By selecting the 



density in the 170 to 200 km height range, the density consists 
primarily of 0 + ions. A positive and approximately direct 
correlation is exhibited by the data, which were obtained during 
solar maximum. 6 

As was discussed earlier, the ionopause height for moderate 
values of the solar wind dynamic pressure depends on a balance 
between the ionospheric particle pressure and the component of the 
solar wind dynamic pressure normal to the ionopause. Since the 
ionospheric particle pressure depends on the product of and the 
sum T e + Tj, we might logically assume that the ionopause height 
will depend on the magnitude of the solar EUV flux in addition to 
the solar wind dynamic pressure. Brace et al. [1990] have recently 
attempted to establish a correlation between the average electron 
density within .several altitude intervals in the central Venus umbra 
and the solar wind dynamic pressure and between the electron 
density and the solar EUV flux. Over a solar cycle the expected 
correlation was obtained; the average nightside electron density 
was positively correlated with the EUV flux and negatively 
correlated with the solar wind dynamic pressure. With the data 
base restricted to solar maximum and low altitude, the nightside 
electron density was negatively correlated with the solar wind 
dynamic pressure, as expected, but showed a weak negative 
correlation with EUV flux which was not expected 
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Fig. 36. Correlation of normalized nightside ion density and normalized 
ionopause height at solar maximum (adapted from [Miller and Knudsen, 
1987). The ion density, measured by the ORPA in the altitude interval 170- 
200 km, is predominantly 0 + and is normalized by dividing it by the 
empirical model density of Theis cl al. [1984]. The measured ionopause 
altitude is normalized by dividing it by the median altitude modeled by 
Knudsen etal. [1982b]. 

THE SOLAR MINIMUM IONOSPHERE 

Strong Solar Cycle Control of the Ionopause Altitude 

The variation of the dayside ionospheric peak density with solar 
cycle phase consistent with Chapman's theory was to be expected 
and has been borne out by measurement [Bauer, 1983; Kliore and 
Mullen, I989J. What was not anticipated or appreciated until 1987 
was the rather dramatic depression of the dayside Cf layer and the 
consequent depletion of the the nightside 0 + layer during solar 
minimum caused by the simultaneous reduction of the average 
solar EUV flux and increase of the average solar wind dynamic 
pressure [Knudsen et al M 1987]. In restrospect, it is rather surprising, 
that this process was not recognized from earlier clues. 

One clue was the discrepancy in nightside density at and above 
200 km altitude between Venera 9 and 10 radio occultation profiles 
obtained during solar minimum and PVO in situ measurements 


made at solar maximum [Savich et al., 1982; Osmolovskii et al., 
1984]. The discrepancy, which exceeds an order of magnitude, is 
illustrated in Figure 37. The densities measured at or near the 
nightside peak by both Venera and PVO radio occultation 
experiments and PVO in situ experiments were comparable. 
Honoring the discrepancy as real, Knudsen et al. [1986] suggested 
that the difference could be understood if at solar minimum the 
transterminator 0 + flux is reduced by an order of magnitude from 
that measured at solar maximum and the electron impact ionization 
produced by suprathermal electrons remained the same. With these 
conditions satisfied, the 0 + layer density would be reduced a factor 
of ten, but the Of layer would be reduced only a small amount. 
However, they failed to recognize how steady shutoff of the Cf 
transport was being effected. 

A final clue came in the form of PVO radio occultation profiles 
in the SZA range 60° - 70° at solar maximum and solar minimum. 
Figure 38 illustrates the average electron density profiles derived 
from several radio occultation profiles measured at solar maximum 
(1980) and at solar minimum (1986). The solar minimum average 
density profile is formed from 10 individual profiles lying between 
60°and 70° SZA. The solar maximum average profile is formed 
from 6 individual profiles falling between 60°and 70° SZA. The 
average solar wind dynamic pressure for the 10 solar minimum 
profiles is 4.8 nPa and is 5.1 nPa for the 6 solar maximum profiles. 
If we assign the ionopause altitude as that at which the density 
decreases below IxlO 4 cm' 3 , the ionopause altitude, as inferred 
from Figure 39, would be at 225 km at solar minimum and at 
575 km at solar maximum. The solar maximum density at and 
above 250 km is more than an order of magnitude larger than that at 
solar minimum. 

Knudsen et al. [1987] attributed this solar cycle change to a 
decrease in the ionospheric pressure by a factor of 3 from solar 

VENUS NIGHTSIDE IONOSPHERE 



Fig. 37. Comparison of the solar maximum median nightside ion density 
altitude profile measured by the PVO ORPA, the solar minimum average 
peak electron density measured by the PVO radio occultation experiment 
[A. J. Kliore, private communication, 1987], and upper limit of the electron 
density at and above 200 km altitude at solar minimum inferred from 
Venera 9 and 10 radio occultation profiles JSavich et al., 1982]. 
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Fig. 38. Average electron density altitude profiles derived from PVO radio 
occultation profiles recorded at solar maximum and solar minimum (from 
[Knudsen et al., 1987]). Each average is derived from several profiles all of 
which were recorded between 60°and 70° SZA. 



Fig. 39. The median ionospheric thermal pressure P, measured at solar 
maximum measured by the PVO ORPA and the inferred thermal pressure at 
solar minimum P/3 (from [Knudsen et al., 1987]). For constant normal 
solar wind pressure, the median ionopause height decreases from 650 km at 
solar maximum to below 250 km at solar minimum. 


maximum to solar minimum. Figure 39 illustrates the median 
ionospheric thermal pressure P, measured by the PVO ORPA at 
solar maximum at 65° SZA. The dashed curve is the predicted 
curve that would be observed at solar minimum in the absence of a 
solar wind. The altitude of pressure balance changes from 650 km 
at solar maximum to 250 km at solar minimum for a constant 
normal component of solar wind dynamic pressure equal to 1 nPa. 
A reduction of the solar maximum pressure profile by a factor of 
three at solar minimum is reasonable because of its nonlinear 
response to the reduction of the solar EUV flux by a factor of two 
or more. The 0 + density, which is directly proportional to the 0 + 
ion production rate, would be decreased a factor of three by a 40% 
decrease in the neutral O density [Keating et ah, 1985] and by a 
50% decrease in the EUV flux. The reduction of the ionospheric 
pressure by a factor of three should, of itself, be sufficient to cause 
the observed reduction of typical ionopause height. However, the 
solar wind also contributes to the decrease. The solar w ind average 
dynamic pressure at solar minimum (1985-1987) was 
approximately 50% larger than it was during solar maximum 
(1979-1981) according to the data shown in Figure 40. 

That the ionopause would be compressed to its limiting altitude 
most of the time at solar minimum may be readily understood by 
referring back to Figure 8. If the peak effective ionospheric 
pressure is reduced by a factor of three, and the center of the solar 
wind dynamic pressure histogram intervals is moved to the right by 
a factor of 1.5, (I am assuming the histogram retains its solar 
maximum form), the percentage of time that the solar wind 
dynamic pressure is less than the ionospheric peak effective 
pressure changes from approximately 78% at solar maximum to 
something less than 10% al solar minimum. 

The Dayside Ionosphere 

No in situ measurements of the solar minimum dayside 
ionosphere exist, but from the available solar minimum radio 
occultation profiles and solar maximum PVO periapsis passes that 



Fig. 40. 100-day averages of the solar wind dynamic pressure pr 2 over a 
solar cycle measured in the vicinity of the earth (adapted from [Lazarus and 
Belcher, 1988]). The dynamic pressure is approximately 50% higher at solar 
minimum (1976, 1986) than at solar maximum (1980). 
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occurred during periods of high solar wind dynamic pressure, we 
can come up with a fairly certain picture. First of all, the dayside 
ionosphere is highly compressed most of the time. No extensive 
diffusive equilibrium 0 + layer overlies the 0 2 + layer such as was 
the case under typical solar maximum conditions. It is reasonable 
to assume that the ion composition is similar to that illustrated in 
Figure 4 for orbit 412, an orbit in which the solar wind dynamic 
pressure was substantially larger than the maximum ionospheric 
particle pressure. Although 0 + ions will be the dominant ion 
constituent above an altitude of about 180 km, the scale height of 
the 0 + ion density will be much smaller than the diffusive 
equilibrium scale height [Mahajan and Mayr, 1989]. 

The Nightside Ionosphere 

With the dayside O* layer effectively removed at solar 
minimum, the nightside 0 + layer should also be effectively 
removed since the nightside 0 + layer is almost exclusively 
produced by 0 + transport. The discrepancy in nightside density at 
200 km altitude between that measured by the Venera 9 and 10 
radio occultation experiments at solar minimum and that measure 
by PVO in situ instruments at solar maximum and illustrated in 
Figure 37 was the first evidence of the missing nightside Of layer. 
Additional evidence is illustrated in Figure 41. Medians of ion 
density derived from ORPA measurements at solar maximum and 
at solar minimum exhibit more than an order of magnitude 
difference in density when extrapolated to a common altitude of 
1800 km. The inference is that a large reduction in the nightside 0 + 


P-V ORPA CONNECTED MEDIANS OF 51 



Fig. 41. Nightside median ion density profiles measured at solar maximum 
and solar minimum by the PVO ORPA experiment (from [Knudsen, 1988)). 
The curves do not overlap because of PVO periapsis altitude changes 
throughout the solar cycle (Figure 3). The ion in parentheses is known or 
thought lo be the most abundant minor ion at the indicated altitude. 


layer density has occurred between solar maximum and solar 
minimum. Brace el al. [1990] have verified that the average plasma 
density above 2000 km altitude at solar minimum is about 10 cm 
or less. 

Brace et al. [1987] have shown that the nightside solar minimum 
ionosphere above about 2000 km altitude is very structured and 
consists of filaments ("tail rays") with a plasma density of a few 
hundred cm' 3 imbedded in a background plasma ("troughs") with a 
density of the order of 10 cm' 3 , as shown in Figures 42 and 43. The 
filaments are thought to be elongated with their long axis more or 
less parallel to the sun-Venus axis. Figure 43 illustrates the 
variation of electron density N c , magnetic field magnitude B and 
component vectors B x , B y and B lt suprathermal O density O sl , 
fast ionospheric ions O f + , electron temperature T e , and total, plasma 
and magnetic pressures P T , P P , and P B through one nightside 
passage recorded in Jan 1984 shortly before solar minimum (see 
Figure 3). Suprathermal ions O s , + appear to compose most of the 
ion density in the rays and troughs and are thought to have a 
temperature in the range 9-16 eV. The fast 0 + ions have an energy 
of approximately 40 eV and comprise of the order of 2% or less of 
the total density [Kasprzak et al., 1987]. The electron temperature 
measured by the PVO electron temperature probe (OETP) is 
approximately 0.5 eV within the rays, a temperature similar to that 
in the solar maximum nightside ionosphere at an altitude of 1000 
km. The electron temperature in the troughs is approximately one 



RU 4432 3545 2*31 2362 2203 237* 2859 3581 4476 

UT 59 4* 35 20 3 -13 -28 -41 -52 DC* 



AIT 4445 3575 2880 2429 2282 24C3 2943 3661 4545 Kn 

Ifk T 57 46 33 IB 1 -H -29 -42 -53 DCC 

Fig. 42. Structure in the solar minimum nightside electron density at high 
altitude (adapted from Brace et al. [1987]). The high density intervals are 
interpreted to be filaments or rays with long dimension approximately 
parallel lo the sun-Venus axis. 
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Fig. 43. Composition and some thermal properties of of the solar minimum 
high altitude tail region(adapted from Brace et al. 119871). The high density 
intervals are interpreted as tail rays. The electron temperature and magnetic 
field strength are low in the rays and high in the troughs. The plasma 
pressure P T and magnetic field pressure P B are in approximate balance at the 
ray edges. The rays and troughs arc believed to be composed primarily oT 
O ions (O**) and are inferred to have a temperature between 9 and 16 eV. 
A component of the 0 + ions O f + has an energy of approximately 40 eV and 
typically composes less than I % of the ions. 


order of magnitude larger. B is strongest in the troughs and is 
parallel to the sun-Venus axis. 

Vaisberg et al. [1976] have reported Venera 10 measurements of 
cool ions at about 1800 km altitude with temperatures of a few eV 
converging toward the sun-Venus axis as shown in Figure 44. The 




Fig. 44. Direction of flow, temperature, and velocity of ions measured by 
the RIEP instrument on Venera 10 in Oct. 1975. Location of the shock and 
boundary of obstacle and flow direction (solid lines in upper portion) are 
from HD model of Spreiter et a!. [1980], Dashed arrows show the low 
energy ion flow direction. The temperature and velocity of the low energy 
component is indicated by dots in the middle and bottom part of the figure, 
respectively (adapted from Russell and Vaisberg [1983]). 

ion temperature is of the order of 1-5 eV and the velocity, if the ions 
are assumed to be 0 + , is 25-50 km s* 1 . The typical number density 
of the heavy component is of the order of 0.5 cm' 3 . Russell and 
Vaisberg [1983] have interpreted these ions as having been 
introduced into the solar wind plasma as it flowed toward the 
terminator within the plasma mantle. It would appear from the 
evidence available at this time that the cool, few eV temperature 
ions observed by Vaisberg et a]. [1976] and Brace et al. [1987] are 
largely 0 + ions produced by photoproduction or charge exchange in 
the mantle layer above the solar minimum ionopause and carried 
into the solar minimum wake. The filaments may have arisen from 
magnetic flux tubes flowing closer to the ionopause than those tubes 
forming the troughs. 

Figure 45 illustrates a nightside solar minimum plasma density 
profile inferred from Venera and PVO radio occupation profiles 
and PVO in situ measurements. The Of peak density is less than 
it is at solar maximum as a result of the elimination of the solar 
maximum 0 + transport ionization source and a possible small 
reduction in the precipitating suprathermal electron flux. Because 
the 0 2 + peak density varies only as the square root of the ionization 
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Fig. 45. Nightsidc median plasma density height profile at solar minimum 
inferred from Venera and PVO radio occultation measurements and PVO 
ORPA in situ measurements (from [Knudsen, 1988]). 


rate, it is relatively insensitive to changes in the ionization rate. The 
solar maximum topside 0 + layer is largely missing and is replaced 
by filaments and troughs of plasma containing 0 + ions and 
presumably H + ions of comparable concentration with ion 
temperatures of the order of 10 eV. The density, composition and 
temperature of the topside plasma suggest that the plasma has 
originated as solar wind plasma flowing within the mantle region 
and subsequently converging into the near-Venus wake. The Ct 
and other ions of Venus origin are added through photoionization, 
charge exchange and possibly diffusion (turbulence) as the solar 
wind flux tubes traverse the dayside mantle region. 

Concluding Remarks 

PVO and previous NASA and Russian experiments have 
revealed a Venusian ionosphere dramatically different than that of 
the earth. The Venusian ionosphere, unprotected and unconstrained 
by a strong internal magnetic field, exhibits a gross morphology and 
many processes and features, such as an ionopause; a supersonic, 
anti sun ward-directed ion wind; and magnetic flux ropes, that are 
substantially different from anything observed in the earth’s 
ionosphere. Our reservoir of knowledge concerning the behavior of 
ionospheric plasma when subjected to constraints different from 
those imposed at earth has greatly increased as a result of these 
successful experiments, and we will be better able to predict the 
behavior of plasmas in new environments. 

Although we now have a large body of experimental data 
defining many of these new processes and features, as is often the 
case, we need additional, refined measurements to discriminate 
between or verify current hypotheses as to the physics behind the 


processes. Some inadequacies in our present measurements needed 
for a better understanding of Venus processes include high spatial 
resolution suprathermal (lOP-lO 5 eV) electron and ion mass 
spectrometer three dimensional flux measurements, and high 
spatial resolution thermal ion bulk velocity measurements. 
Hopefully, PVO and the Venera missions will not be the last plasma 
physics missions to Venus. 
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Abstract . Suprathermal electron density, tem¬ 
perature, integral flux above 58 eV energy and the 
magnetic field are presented for two Pioneer-Venus 
orbits passing through the antisolar region of the 
Venus ionosphere. One of the orbits passed 
through an ionospheric hole. The second orbit did 
not exhibit a hole signature and is representative 
of other orbits examined in the antisolar region. 
The suprathermal electron properties exhibited by 
both orbits are consistent with median values and 
their variation reported for the Venus umbra in a 
previous study. The suprathermal electron temper¬ 
ature, density, and integral flux within the hole 
region are similar in magnitude to these same 
quantities observed exterior to the hole and give 
no suggestion of the existence of a significant 
parallel electric field. We infer that iono¬ 
spheric holes are formed by a less energetic pro¬ 
cess than the parallel electric fields responsible 
for aurorae in the auroral zone of the earth. We 
favor the more benign process of plasma loss 
through chemical recombination at low altitude as 
dayside ionospheric plasma with imbedded, weak 
magnetic field flows toward the antisolar region. 
The magnetic field intensifies as the plasma den¬ 
sity within a tube decreases to maintain pressure 
balance. The suprathermal electron energy flux 
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entering the atmosphere within the hole system is 
estimated at 2.4 x 10 3 erg cm~^ s - ^ and will 
produce an ionospheric layer with peak density of 
the order of lxl0 4 cm~ 3 at an altitude of approxi¬ 
mately 145 km. The layer, composed almost entire¬ 
ly of molecular ions, will decrease in density 
rapidly with altitude, a property exhibited by 
holes. 
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1. Introduction 

An unusual ionospheric phenomenon observed in 
the Venus nlghtside ionosphere has been termed a 
"hole." A hole is a region in which the electron 
concentration is strongly depleted and the mag¬ 
netic field is increased in strength and is ap¬ 
proximately radial in direction (Brace et al., 
1982, Luhmann et al., 1982). The purpose of this 
paper is to present suprathermal electron proper¬ 
ties measured in a hole by the Pioneer—Venus re¬ 
tarding potential analyzer (RPA), properties not 
previously published, and to use this new infor¬ 
mation to limit the range of proposed mechanisms 
of hole formation. 

Grebowsky and Curtis (1981) suggested that 
holes may be a manifestation of parallel electric 
fields similar to those observed in the terres¬ 
trial auroral ionosphere. More recently, Grebow¬ 
sky et al. (1983) have developed a simple two- 
dimensional model in which holes develop as plas¬ 
ma, transported from the dayside ionosphere, is 
lost into the lower ionosphere through normal 
chemical processes. The magnetic field within the 
tube concurrently intensifies to maintain pres¬ 
sure balance. Elements of this process were 
previously suggested by Knudsen et al. (1981, 
1982) and Brace et al. (1982). 
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We present measurements of suprathermal elec¬ 
tron properties along two periapsis passes of the 
Pioneer-Venus spacecraft through the antisolar 
region of Venus. The spacecraft passed through 
an ionospheric hole on the one pass and did not 
encounter a hole on the second. The suprathermal 
electron properties measured through the hole re¬ 
gion do not indicate the presence of a parallel 
electric field and are typical of properties 
measured through non-hole regions. We infer from 
the results of this study and previous studies 
that ionospheric holes can be formed by processes 
other than parallel electric fields and that it 
is unlikely that significant parallel electric 
fields play a role in forming ionospheric holes. 

2. Method of Measurement 

The suprathermal electron quantities which we 
report herein were measured by the Pioneer-Venus 
retarding potential analyzer while operating in 
its Buprathermal electron mode (Knudsen et al., 
1979, 1980). In this mode the instrument records 
the integral electron current as the retarding 
potential sweeps from 0 to -50 V. The instrument 
measures the background current at -58V. Since 
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the vehicle potential V gc is typically of the 
order of —1.5 V negative relative to the plasma 
when the spacecraft is within the ionosphere and 
in the umbra, the background current is measured 
at an effective retarding potential of approxi¬ 
mately -60 volts. We have shown in a recent 
study that the electron current-voltage (i-V) 
curves may be approximated quite accurately by 
.assuming the electron gas to be composed of two 
Maxwellian components with concentrations Nj and 
N 2 and temperatures Tj and T 2 (Knudsen et al., 
1984). and Tj approximate the thermal elec¬ 
tron density and temperature, respectively. N 2 
and T 2 are the density and temperature of a 
suprathermal electron component and typically have 
values equal to 3 cm -3 and 14 eV, respectively 
(Knudsen et al., 1984; Gringauz et al., 1979). 
The suprathermal electron distribution has an 
energetic tail which may be represented by a third 
Maxwellian distribution with a typical density Ng 
of 'x-O.l cm J and a temperature of 100 eV (Gringauz 
et al., 1979; Knudsen et al., 1984). This more 
®^®tgetic component is primarily responsible for 
the background integral current per unit area i co 

JO 

measured with a retarding potential of -58 V. We 
shall indicate the variation of N 3 by presenting 
the variation of i^g normalized by tt and e, the 
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electronic charge. i 5g /(Tr e) is an approximation 
of the integral electron flux per steradian above 
60 eV energy. The exact expression is (i 5g 
“ V di/dV)/(iTe). Both terms in the expression are 
comparable in magnitude. 

A more detailed discussion of the least- 
squares method we use in deriving the parameters 
N 2 S Tj, T 2 and N g from an RPA i-V curve and an 
example of an i-V curve with least-squares fit has 
been published previously (Knudsen, et al., 1984). 
We summarize here some of the uncertainties dis¬ 
cussed in more detail in the previous paper. The 
RPA electrometer root-mean-square current noise is 
6 x 10 A and yields an uncertainty in the 
quantity i^/Cire) of approximately 1 x 10^ cm"^s" 1 
We have estimated that the uncertainty in N 2 
and T 2 caused by random noise in the RPA is ap¬ 
proximately 25% for the median values of these 
quantities, N 2 - 3 cm -3 and T 2 - 14 eV, observed 
in the Venus umbra. The uncertainty is estimated 
to Increase to 100% for N 2 equal to or less than 
approximately 0.5 cm -3 . 

The vehicle potential within the holes and in 
other regions where the thermal plasma (TM 
eV) density is very small (<_10 2 cm" 3 ) is un¬ 
known. It may become of the order of -T 2 /e where 
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T 2 is expressed in eV. If the vehicle potential 
does become this negative in these regions, the 
value of N 2 derived with the assumption that it 
is a uniform -1.5 volts will be too small by 
a factor of exp(-T 2 /1.5 eV) where T 2 is in eV. 

An occasional N 2 - T 2 pair may exceed the 25% 
uncertainty estimated above because the i - V 
curve from which it is derived by least—squares 
fitting has one or more erroneous current values 
caused by telemetry dropout or an electrometer 
switching transient. These bad current values 
can be recognized by inspection but the data 
reduction software is not sufficiently sophisti¬ 
cated to recognize all of them. 

3. Experimental Observations 
Orbit with Hole 

In Figure 1, we have plotted the thermal elec¬ 
tron concentration measured by the P-V Langmuir 
probe N g (Krehbiel et al., 1980), the supra- 
thermal electron concentration N 2 with tempera¬ 
ture T 2 and the integral current i^g. In Figure 
2, the total magnetic field |b| and the com¬ 
ponents B x , By, and B z measured by the P-V magne¬ 
tometer (Russell et al., 1980) are presented for 
the same time period. The coordinate system is 
centered in Venus with the z axis perpendicular 
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to the earth s ecliptic plane and directed 
northward. The sun is contained in the x-z plane 
with positive x toward the sun. At approximately 
1:39 UT, the thermal electron density N e is 
strongly depleted relative to the surrounding 
plasma density, and the magnetic field is en¬ 
hanced and is primarily in the +B X direction 
(toward the sun). These two signatures are those 
characterizing an ionospheric hole. We observe 
that i^g /( ne) does not reflect any significant 
enhancement as the spacecraft passes through the 
hole region. N 2 is somewhat depressed relative 
to the surrounding values. This decrease within 
the hole may be only apparent because, as 
explained above, the spacecraft potential may be 
more negative within the hole than in the 
surrounding ionospheric region. Two values of T 2 
appear somewhat enhanced - 20 eV instead of the 
typical 14 eV — on the walls of hole region and 
have correspondingly small values of N 2 « 
isg/^e) reveals no increase in flux above 60 eV 
energy. At the present level of our investiga¬ 
tion, we are unable to attach any significance to 
the occurrence of these two values. The least- 
squares analysis software yielded a standard 
deviation of 4 eV for the two larger tempera- 
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tures and 1 eV for the remainder which indicates 
the two i-V curves were somewhat noisier than the 
remainder. We interpret these measurements of 
i 58‘ n 2» an< * t 2 as indicating no significant 
enhancement of electron flux or energization of 
electrons within the ionospheric hole region. 

In Figure 3 we have expanded the time interval 
near UT 1:39 and plotted the cosine of the av¬ 
erage pitch angle 0 between the RPA normal (posi¬ 
tive in the outward direction) and B. The direc¬ 
tion of "F has been calculated only at UT 1:38:30 
and is assumed to be approximately constant over 
the interval 1:37 to 1:40 UT. We have plotted 
the thermal electron density measured by the 
RPA instead of N g measured by the Langmuir probe 
because of the exact time coincidence of the N 1# 
N 2 , T 2 , and cos(0 ) values of the RPA. Nj is 
calculated with the assumption that the vehicle 
potential is a constant -1.5 V relative to the 

plasma potential. N, differs in value from N 

x e 

because it has not been corrected for the factor 
exp(-eV sc /kT 1 ), but it serves to locate the 
strong thermal depletion region. The objective of 
Figure 3 is to reveal any pitch angle anisotropy 
of the electron flux that may be detectable be¬ 
cause of the changing orientation of the RPA with 
spacecraft spin. Although the P-V RPA is a 



12 


planar instrument, its "effective” acceptance 
solid angle is approximately ir steradians, a 
cone with half-angle of approximately 60°. The 
RPA normal describes a cone with half-angle of 
25° and with axis of the cone directed toward 
the south ecliptic pole as the spacecraft ro¬ 
tates. Since S’ is approximately parallel to the 
earth ecliptic plane during the time interval of 
interest, the RPA primarily samples the electron 
flux with pitch angles nearer to 90° than to 0° 
and 180°. The RPA is programmed to record an i-V 
curve at four celestial longitudes separated by 
approximately 90°. This cyclical sampling pat¬ 
tern is reflected in the cos(0) curve of Figure 3 
with the exception that the instrument did not 
record (or transmit) the i-V curve at approxi¬ 
mately 1:37:15 UT, the time at which cos(6 ) would 
have been negative. 

The cyclical variation of i 5 g/( tt e) between 
1:37:30 and 1:40:30 UT suggests that some elec¬ 
tron flux anisotropy exists in the flux with 
energy greater than 58 eV. When cos( 0 ) is 
negative, i 58 /( Tfe) exhibits a maxima, and when 
cos(0 ) is most positive, i 5 g/(*e) exhibits a 
minima. Although the uncertainty in i 5 g/(u e ) 
resulting from the electrometer noise level is 
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comparable in amplitude to the cyclical varia¬ 
tion, we believe the cyclical variation to be 
indicative of a somewhat greater electron flux 
with energy above 58 eV traveling down the magne¬ 
tic flux tube toward the atmosphere than that 
traveling up the flux tube. 

There is no clear indication in Nj or N 2 of 
flux anisotropy in these lower energy electrons. 
Tj is approximately 1 eV within the hole region. 
These variables show little variation through the 
hole region except for the two larger values of 
T 2 and the corresponding small values of N 2 at 
approximately 1:38 and 1:30 UT. 

Orbit without Hole 

In Figure 4 we illustrate the thermal electron 
density from the Langmuir probe and suprathermal 
electron properties from the RPA for orbit 
1197. Figure 5 illustrates the magnetic field 
for the orbit. The thermal electron density in 
orbit 1197 is more irregular and depressed below 
that of orbit 1203 and is essentially absent from 
UT 1:38 to 1:43. This latter region is not 
interpreted as an ionospheric hole because the 
field is not enhanced over the general level of 
the field elsewhere and is not restricted to 
being largely in the radial (B x > direction. We 
interpret the magnetic field behavior as that of 
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the interplanetary magnetic field draped over the 
front of Venus and extending into the magnetic 
wake with the field in opposite directions on 
either side of a neutral sheet. The general 
level of the background electron current measured 
at -58 V retarding potential is similar to that 
of orbit 1203. Also, i^g/(ire) does not exhibit 
any unusual behavior in the UT interval, 1:38 to 
1:43, in which N g is strongly depleted. The 
suprathermal electron temperature T 2 exhibits 
some large values 050 eV) in regions where N 2 is 
of the order of 1 cm These may be noise 

values resulting from the least squares fitting 
process. The background flux i^g/(^ e) does not 
reflect any unusual corresponding variation. 

We have estimated the total energy flux de¬ 
posited in the atmosphere by the N 2 and Ng Max¬ 
wellian populations by assuming that the flux is 
isotropic. The energy flux passing through unit 
area from one side in one second for a Maxwellian 
distribution is 

1) Energy flux ■ (2/CTTm)) 1 ^ 2 N (kT)"^ 2 
where the symbols have their usual meaning. Us¬ 
ing the measured values of Nj, Tj, N 2 , T 2 , the 
largest energy flux we computed for the UT time 
interval 1:37:30 to 1:39:40 is 2.4xl0“ 2 erg cm -2 
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sec *. The largest energy flux In the depressed 
N g interval of orbit 1197 is 3.4xl0 -2 erg cm -2 
sec - *. 


4. Discussion and Conclusions 

We interpret the experimental observations pre¬ 
sented in Figures 1-5 as evidence that an ener¬ 
getic process such as a significant parallel 
electric field formed along the radial magnetic 
field of an ionospheric hole is not responsible 
for the hole formation. As the Pioneer-Venus 
spacecraft moved through the hole region in orbit 
1203, the integral current i^g measured with 
retarding potential at -58 V reflected no in¬ 
crease in flux. At earth, a spacecraft equipped 
with a low energy electron detector, such as the 
soft particle spectrometer on the ISIS 2 space¬ 
craft, shows a definite and substantial increase 
in total electron flux, energy flux, and mean 
energy of the flux as the spacecraft passes 
through an auroral feature (Winningham and Heik- 
kila, 1974; Knudsen et al., 1977). A peak in the 
differential energy flux develops within the 
auroral precipitation region at an energy of a 
few KeV as a result of the electrons dropping 
through the parallel electric field formed at 
higher altitude. Although the P-V RPA is incap- 
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able of measuring the differential electron flux 
at energies above approximately 45 eV, it would 

reveal any significant increase in integral flux 
above 58 eV. 

The differential flux below 50 eV is also in¬ 
creased within an auroral form as a result of the 
production of secondary electrons by the KeV pri¬ 
maries. We observe no increase in N £ which would 
signal an increase in the production rate of sec¬ 
ondary electrons. 

The energy flux of the suprathermal electrons 
within the ionospheric hole of orbit 1203 has 
been estimated above at 2.4xl0 -2 erg cm -2 sec -1 . 
This energy flux is one to two orders of magni¬ 
tude less than that producing faint, diffuse 
aurora (IBC1) on earth (Chamberlain, 1961). It 
is also an order of magnitude less than the 
typical electron energy flux entering the Earth's 
cusp region (Heikkila and Winningham, 1971). 

The presence of a suprathermal electron flux 
within hole magnetic flux tubes with magnitude 
typical of that observed over the nightside Venus 
ionosphere contributes to our understanding of an 
observed property of holes. When the Pioneer- 
Venus. spacecraft passes through a hole on both 
the inbound and outbound legs of a single periap- 
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sis pass, It passes through the hole on the 
inbound leg at a lower altitude than the hole on 
the outbound leg. The plasma density is observed 
to be less depleted in the inbound hole, which is 
at lower altitude, than in the outbound hole 
(Luhmann et al.,. 1982). The explanation which 
has been given previously for this observation is 
that collisions of ions and electrons with neut¬ 
rals at the lower altitudes permit diffusion of 
plasma into the hole region from the surrounding 
denser plasma. We point out here that the supra- 
thermal flux observed within the hole of orbit 
1203 is sufficient to maintain a layer of molecu- 

lar ions with peak density of the order of 1x10^ 
-3 

cm at an altitude of approximately 145 km 
(Spenner et al., 1981; Cravens et al., 1983). 
The 0 + density produced by the electron flux is 
less than approximately 1x10^ cm“^. Hence, if we 
assume that the 0 + ions initially contained in a 
flux tube have been lost largely through downward 
diffusion and recombination, a plasma layer with 
a topside molecular ion Beale height will still 
be maintained with near-typical peak density at 
an altitude of approximately 145 km by the supra— 
thermal electron flux. 
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We have argued that the suprathermal electron 
properties measured In one clear example of an 
ionospheric hole do not reflect the existence of a 
significant parallel electric field accelerating 
electrons toward Venus. We have shown that the 
suprathermal electron properties within the hole 
region and in the adjacent non-hole region are 
similar to suprathermal properties measured in an 
orbit passing through the same general antisolar 
region six days previously. The properties ex¬ 
hibited in both orbits are consistent with median 
properties and their variation derived in a pre¬ 
vious study of the low altitude Venus umbra (Knud- 
sen and Miller, 1984). The present study, the pre¬ 
vious study, and a cursory examination of the 
suprathermal electron properties in several orbits 
around 1197 and 1203 all suggest to us that the 
Venus nightside ionospheric holes are not formed 
by parallel electric fields of any significant 
magnitude. We are thus led to favor a less ener¬ 
getic process for forming the holes. 

Insufficient data are available to be very 
specific about the process or processes forming 
the holes. We do not know the two-dimensional 
shape of the holes — let alone the three-dimen¬ 
sional shape (Brace et al., 1982; Marubashi et 
al., 1985). Neither do we know the range of envi- 
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ronmental conditions under which the holes form or 
do not form. We do not attempt to be very spe¬ 
cific, therefore, but one less energetic process 
which we suggest is involved in the formation of 
holes and to which several authors have con¬ 
tributed (Knudsen et al., 1981, 1982; Brace et 
al., 1982, Grebowsky et al., 1983) is as follows: 
-As high beta plasma from the dayside ionosphere 
flows supersonically toward the antisolar region, 
a portion of the plasma is continuously lost 
through a downward component of flow into the low 
altitude chemical sink where atomic ions are con¬ 
verted to molecular ions which then recombine. As 
the plasma density decreases within a magnetic 
tube, the magnetic field increases to maintain 
pressure balance across the tube wall. A hole 
forms when the loss process is sufficiently rapid 
compared with the horizontal transport time. 
Slowing of the horizontal plasma flow by passage 
through a shock (Knudsen et al., 1980) may be 
essential in forming the holes in which the den¬ 
sity is most strongly depleted (Grebowsky et al., 
1983). 
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Figure Captions 

Figure 1. Thermal and suprathermal electron properties measured as 
the P-V spacecraft passed through an ionospheric hole on the 
nightside of Venus. N g is the thermal electron density and N 2 
and T 2 the suprathermal electron density and temperature, 
respectively. i 5g is the integral current measured by the P-V 
RPA with retarding potential at -58 V. 

Figure .2. Magnetic field measured along the same orbit segment as that 
in Figure 1. The z coordinate axis points toward the North eclip¬ 
tic pole and the x-z plane contains the Venus sun vector with 

positive x toward the sun. 

Figure 3. Thermal and suprathermal electron properties measured in 
the Immediate vicinity of the ionospheric hole of orbit 1203. 
is an approximate thermal electron density measured by the P-V 
ORPA; is the angle between the RPA outward pointing axis and 
the magnetic field vector. The integral current i 58 variation 

suggests the electron flux with energy above 58 eV may be 

anisotropic with the greater flux directed toward the Venus 
atmosphere• 

Figure 4. Thermal and suprathermal electron properties measured as 
the P-V spacecraft passed through a region of the nightside 
ionosphere not exhibiting an ionospheric hole. 

Figure 5. Magnetic field measured along the same orbit segment as 
that in Figure 4. 
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I. BACKGROUND 

The Unified Abstract Data System, UADS, was developed as a computer based, 
interactive system to provide Pioneer Venus Orbiter (PVO) experimenters quick 
and convenient access to data being acquired by the mission. Data were 
entered into a central facility by each experimenter through local terminals; 
the UADS software then sorted and cross-filed the data by orbit numbers. Eacn 
experimenter could then access any information from the central facility 
through his local terminal. The UADS proved invaluable in providing quick 
access to cross-compare results between experiments, especially in the early 
phases of the mission when rapid feed-back to Mission Control was critical. 

The UADS further produced data files on magnetic tape, suitable for 
archival purposes with the National Space Science Data Center, NSSDC. The 
files so produced contained Low Frequency Data (LFD), Special Event Data 
(SED), and Composite Data (CD). Both the SED and CD record formats were 
defined and controlled by each experimenter. The LFD combined experimental 
results and orbital information in an unchanging, predefined format. The 
experimenters were responsible for entering data peculiar to their instrument, 
while the Pioneer Mission Office (PMO) supplied orbital parameters and 
experiment pointing information through a subset of the Supplementary 
Experimenter Data Record (SEDR). 

The software for the UADS was developed, under contract, for this purpose 
by Bendix Corporation and the computing resources of Tymshare Corporation were 
purchased to serve as the central facility. 

As the mission progressed, the need for rapid access to data, cross- 
referenced by orbit, decreased, and under the pressing need to reduce costs a 
decision was made by the PMO to terminate the contract with Tymshare Corpor¬ 
ation as of January 1981. Without support from Tymshare Corporation the UADS 
would not be available for preparing NSSDC submissions, so the PMO agreed to 
develop a substitute system to provide archival data records. 

II. PIONEER VENUS ARCHIVAL DATA AT THE NSSDC . ... 

In conformance with NASA policy. Pioneer Venus data are archived at the 
NSSDC, whose charter it is to acquire, store and disseminate space science 
data for further analysis beyond that provided by the principal investigators. 

Pioneer Venus data have been archived with the NSSDC in the format 
produced by the software of the UADS. Various plans have been suggested to 
implement the continued submission to the archives now that the UADS has been 
terminated. Several plans were submitted for consideration to the May 6, 1983 
meeting of the OMOP Committee. These were, in brief, 

1. The experimenters send their data to the PMO, who periodically 
contract with a time sharing vendor to run, as a batch job, the archival 
feature of the UADS. The primary advantage of this would be to produce 
archival data in the same format as that already at the NSSDC with very little 
reconfiguration of any existing system. An important drawback is an estimated 
very high cost (about $30K for each six months of data processing). 




2. The experimenters send their data to the PMO who will arrange to use 
an existing NASA Ames IBM 4341 computer to run a modified version of the 
archival feature of the UADS. The primary advantage of this approach would be 
that it will be relatively inexpensive to use. 

Both plans (1) and (2) retain the tape formats used by the UADS. As is 
discussed below, this was argued to be an important drawback. On the basis of 
this a third plan was suggested, which is as follows: 

3. The experimenters individually send their data directly to the NSSDC. 
Each experimenter would need to create a self-contained submission which would 
Include experiment data, orbital parameters and experiment pointing informa¬ 
tion. This concept would result in NSSDC archives organized by experiment 
rather than by orbit. An advantage of this is that the PMO will no longer 
need to pay for the elaborate and expensive UADS. Further, data requests to 
the NSSDC are usually made for a single instrument; this approach would 
eliminate the need for a requestor to strip out the data of interest from 
records of all instruments. 

This concept would also allow data to be provided to the NSSDC in a tape 
format which is more convenient to use than that provided by the UADS. The 
LFD file is recorded on tape both in 32-bit binary format and in IBM 360 
floating point format; formats which may present conversion problems for a 
user employing a small system. Such a user, especially in a remote location, 
would have great difficulty reading such a tape, or be forced to develop 
conversion software. 

III. COMMITTEE RECOMMENDATIONS 


A committee was appointed to consider and recommend a procedure to 
Implenent the continued submissions of PVO data to the NSSDC. The committee 
consisted of: 


Peter Ford, 

Larry Brace, 
Robert Theis, 
Winifred Cameron 
Pat Barclay, 
and Roger Craig 


MIT 

GSFC 

GSFC 

NSSDC 

Bendix Corporation 
ARC, PMO 


The committee discussed the above and agreed that the production of 
archival files which are easily accessible and usable to a variety of users 
was of high priority. It was agreed that to achieve this it is required that 
files be generated in a style more consistent with the usual NSSDC submis¬ 
sions, namely organized by instrument rather than by orbit. It was further 
agreed that in order to broaden the potential user base the UADS format 
(32-bit binary and IBM 360 floating point) be replaced with a more generally 
used format. 


The results of committee discussions of the above are guidelines for 
submitting PVO data for the NSSDC. The guidelines are designed to: 





1. standardize and improve the data formats for each experiment so that 
the tapes can be easily read by a requestor. 

2. store the data by experiment, 

3. provide for orderly updating of data, 

4. provide for ultimate replacement of existing UADS generated archival 
tapes. 

A new format was designed which is to be used by each experimenter for PVO 
NSSDC data experimenter, and can be easily read by computers expected to be 
accessible to a requester. One major difference between this and the existing 
LFO format is the use of text (ASCII) data formats, eliminating both binary 
and IBM floating point formats. 

The proposed new submission format is self-defining in the sense that the 
first three records on each tape define the data parameters, value representa¬ 
tions, and missing data (file) indicators. The first tape record defines the 
order in which the variables appear in the subsequent data records, in a 
manner similar to that used for SEDR trajectory data (SEDR file 5). The 
second tape record will contain a FORTRAN-compatiole format list describing 
the field sizes and representations of each data value in the order defined in 
record 1. This format may be used to decode all subsequent records on the 
tape. 

The third tape record will define a unique value associated with filler 
(missing) data for all variable fields. It is formatted according to the 
format used in record 2, and is immediately followed by the start of actual 
data records (records 4 and beyond). 

The following is the proposed new format, with examples as applied to OETP 
instructions. 


PROPOSED PIONEER VENUS NSSDC LOW-FREQUENCY DATA FORMAT 


This document describes a suggested format to be used by all investigators for 
the submission of their data to the National Space Sciences Data Center. The 
overall specification will require that all data be coded into ASCII, and 
written onto standard 1/2 inch 1600-bpi 9-track tapes. The logical record 
length will be fixed for a given tape, as well as the physical blocksize. 
Blocksizes should be large enough to avoid wasting tape, but should not exceed 
8000 bytes in order to avoid making excessive demands on user programs for 
memory. The first three records of any of these tapes will be formatted as 
follows: 


'3 



Record 1: The format to be used is (13,n(lX,A4)) where "n" Is the number of 
data items in each record. 


* 

4 

ELTE 

ELNE 

HI 

VS 



(for OETP) 

7 

ETEM 

SPOT 

TONE 

TTWO 

XVEL 

YVEL 

ZVEL (for ORPA) 

T 

T 

T 

T 

T 

T 

T 

t 

3 

5 

10 

15 

20 

25 

30 

35 


Example 1: The first record in each tape file. Note that new value types 
with new 4-character designations can be added as necessary. The date, 
time, orbit and time-tag items are not included in the list, because they 
are common to all data records. 

* number of data items "n" 


Record 2: This record contains the format in which all succeeding records are 
written. The first 4 format items specify the date, time, orbit, and 
time-tag, and will appear in the same format on all tapes. 


(18,19,15,IS,4F9.2) (Appropriate for OETP) 

t 

l 


Example 2: The second record in each tape file 


Record 3: This record will contain zeroes for the first four fields (date, 
time, orbit, and time-tag), and in addition will have a fill value in each 
data value location. This value will be used by any program reading the data 
to identify fill data in subsequent input records. 


0 

0 

0 

0999999.99999999.99999999.99999999.99 

T 

T 

T 

I T t t t 

8 

17 

22 

28 37 46 65 64 


Example 3: The third record in each tape file. (Appropriate for OETP). 






Record 4 to : These records contain the date, time, orbit, and time-tag for 
each time which has any non-fill data. 


1981207 

43527786 

8/9 -1788 

2345.67 

78543.89999999.99 

16.20 


T 

T 

T 

T 

1 

T T 

T 


8 

17 

22 

28 

37 

46 55 

64 



Example 4: All records after the third in a tape file. (Appropriate for 
OETP). 


As can be inferred from the above example, the date is coded as YEAR, DAY OF 
YEAR (1-366) with 19 included in the year. The time is in milliseconds of the 
day, orbit number is self-explanatory, and the time tag is the usual value 
ranging of the day from -1800 to 1800 in increments of 12. 

The project-provided tape of SEDR information would be the source of the 
official dates and times to be used by all other investigators. 

Nothing In the above format would preclude investigators from producing a tape 
containing the data from more than one experiment. 

The external label on the tape should be type-written, and contain the 
following information: 

o Full name of experiment data contained on tape.^ 
o Start date, time, and orbit number of data on the tape. 

o Stop date, time, and orbit number of data on the tape. 

o Production date of the tape. 

o The density (1600-bpi) and number of tracks (9) at which the tape was 
recorded. 

o An estimate of the amount of tape used. 

o The physical blccksize used in writing the tape. 

o A name and phone number of the individual responsible for the tape. 


* Example: "Pioneer Venus Orbiter Electron Temperature Probe". 
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Each Experimenter will be responsible to modify his software accordingly 
and when appropriate, submit data tapes to the NSSDC. As data are updated, 
new tapes should be sent directly to the NSSDC with instructions as to which 
tapes it replaces, or appends. 

Documentation is important in these submissions. All submissions must 
conform to the NSSDC instructions "Guidelines for Submitting Data to the 
National Space Science Data Center." These instructions delineate the kinds 
of documentation which must accompany shipments. 


IV TIMETABLE 

1. Submissions should oe made to the NSSDC for each experiment by 
October 1, 1983 to complete the data set for orbits 1 to 1000. 

2. Submissions should be made to the NSSDC for each experiment by 
November 1, 1983 to replace all data currently on file in the UADS format. As 
replacements are received, the NSSDC will destroy the data in the UADS format. 
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INTRODUCTION 


This document describes the Pioneer Venus (PV) Retarding Potential Analyzer (RPA) National Space 
Science Data Center (NSSDC) low frequency data (LFD) experimenter data record (EDR) tapes. These 
tapes contain one file for each orbit for which data are available. Each file consists of three information 
records followed by n data records. The value of n is variable for each orbit (file), but does not exceed 
721. Each data record consists of four quantities associated with a time tag followed by 25 physical 
quantities unique to the RPA. The formats for the three information records and the data records are 
described in the attached report "INSTRUCTIONS FOR DATA SUBMISSIONS TO THE NATIONAL 
SPACE SCIENCE DATA CENTER, A COMMITTEE REPORT" by Roger A. Craig and dated August 11, 
1983. Ephemeris data for the time tags in each record is contained on a supplementary experimenter data 
record (SEDR) tape supplied to NSSDC by the Pioneer Venus Project Office and will not be described in 
this document. 

The PV RPA is described in considerable detail by Knudsen et al. [1979,1980], The principles of 
measurement are also described therein together with some of the factors affecting the accuracy of the 
derived quantities. Additional information on the theory of measurement by an RPA is presented by 
Knudsen [1966]. 

We describe in this documentation the list of 25 quantities provided in the RPA EDR tapes and their 
limitations appropriate to the submission beginning Nov. 4, 1993. 

The PV NSSDC LFD SEDR tapes have time tags at 12 second intervals from 30 minutes prior to periapsis 
to 30 minutes after periapsis. These time tags are the tags specified in the first four quantities of each of 
the EDR data records. All PV instruments are to report their data at these common time tags for the 
purpose of easy inter comparison of data. Principal Investigators (Pis) with instruments with a sampling 
period much less that 12 seconds are to report the average of measured quantities over a 12 second interval 
centered on the time tags. The RPA, because of a low telemetry word assignment, records at most one 
current-voltage (I-V) characteristic curve per spacecraft spin period. Except for one set of 14 orbits, the 
spin period of the PV spacecraft has been about 12 seconds. Thus, RPA physical quantities are derived at 
intervals of 12 seconds or more. Since the RPA operates in several modes, a particular quantity such as 
thermal electron temperature may be typically measured at much longer intervals. The thermal electron 
temperature is typically measured at either approximately 48 or 60 second intervals. In a few orbits, it was 
measured at 12 second intervals. 

Since the RPA measures its quantities at time intervals of 12 seconds or greater, averaging RPA quantities 
is meaningless. Consequently, RPA quantities are reported in the EDR data records as derived. The 
derived value is, for the most part, reported in that data record whose time tag is within plus or minus 6 
second of the actual time of measurement. (An individual RPA I-V curve is measured in a time interval of 
approximately 0.2 seconds. The average time interval between measurement of successive I-V curves is 
approximately 12 seconds.) The actual time at which the quantity was recorded is given in the record 
immediately following the four quantities specifying the time tag. Thus, an investigator can use his own 
scheme of averaging to assign a value to time tags or may plot or otherwise use the quantities at their times 
of actual recording. When no RPA quantity is available for recording at a specified time tag, no data 
record for that time tag is written on the NSSDC EDR tape. 

RPA quantities may be unavailable for assigning to a specific time tag for several reasons as follows: The 
spacecraft data format in use at the time may not have contained any words for the RPA. The RPA may 
have been turned off for power conservation reasons. The spacecraft telemetry bit rate and/or data format 
may have been such that an RPA I-V curve was recorded only at long time intervals. 
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RPA MEASURED QUANTITIES 


The PV RPA is described together with some of the principles of measurement m some detail by Knudsen 
et al H979 19801 Many of the factors affecting accuracy are also described therem. We present in this 
section the quantities recorded on the NSSDC EDR data files following the four time tag quantities, their 
nominal uncertainty and measurement noise level, and additional limitations of the quantities. 


Table 1 lists the symbol, quantity, measurement range with units in which the quantities are quoted, noise 
level of measurement, and uncertainty of the measurement for the quantities reported by the RPA. 


UT: UT is the universal time in milliseconds assigned to the physical quantities recorded in this record. 
UT will typically, but not always, lie within plus or minus 6 seconds of the time of day assigned to the time 
tag of this record. UT should be accurate to within plus or minus 0.1 second. 

TOT1: TOTI is the total ion density of the plasma in cm-3 and is derived from the FORTRAN expression 


TOTI=FlI/(VN*e*Area) 


where FI I is the first ion current measured with zero retarding potential, VN is the component of the ion 
bulk velocity relative to a coordinate system fixed in the spacecraft and parallel to the RPA axis, e is the 
electronic charge, and AREA is the effective area of the RPA collector. When a lst-squares analysis is no 
possible, VN is assumed to be the component of the spacecraft velocity in ecliptic coordinates parallel to 

the RPA axis. 

H+- H+ is the hydrogen ion density. When the RPA is operating in one of its peaks mode, H+ will be 
detected and recorded only when its density is greater than approximately 10% of the sum of more massive 
ion densities H+ can be the second most abundant ion and still not be recorded when the RPA is operating 
in its two peaks mode. The uncertainty of the H+ density also depends on its density relative to that of 
more massive ions. For an H+ density comparable to that of more massive ions, the accuracy should be of 
the order of 10%. The detection noise level for H+ is estimated at 50 cm-3. Additional discussion of the 
RPA ion peak detection capability and limitation is given by Miller et al. [1984], 

0+- 0+ is the oxygen ion density. It will be detected in the presence of more massive ions only when its 
density is greater than approximately 10% of the sum of more massive ions. The RPA does not resolve 
C+, N+, or 0+. We have assumed in our least-squares fitting that ([C+] + [N+])/[0+] is constant at 0.07, a 
value derived from PV ion mass spectrometer results. 

02+: 02+ is the symbol assigned to the ion with atomic mass of 32. In fitting the RPA I - V curves, it has 
been assumed that an ion of mass 29 is present in the plasma and has a density that is a fixed fraction of the 
02+ ion density at a given altitude. The density ratio of 02+ to the ion of mass 2 varies om 
approximately 5% at 150 km altitude to 90% at and above 200 km. 


C02+: C02+ is the density of the carbon dioxide ion with atomic mass of 44. 

TI: TI is the ion temperature and is assumed to be the same for all ion masses. It is one of the adjustable 
variables in the least-squares analysis of ion sweeps. 

VxVyVz: This vector, specified in the earth ecliptic coordinate system, provides the magnitude 
(MAG =sqrt(Vx**2 + Vy**2 +Vz**2)) of the ion bulk velocity in the direction defined by the unit vector 
(Vx/MAG Vy/MAG, Vz/MAG), and is measured in each ion sweep for which a successful least-squares fit 
has been obtained. An approximate value of the vector ion bulk velocity can be obtained from three 
successive vectors obtained in three successive sweeps since the three unit vectors are not parallel to eac 

other. 
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FI I: FI I is the saturation (first) current measured in an ion I-V sweep. The retarding potential is 
programmed to be slightly negative of plasma potential during this measurement. FI I is measured relative 
to the ion current measured with the retarding potential equal to 37V positive [Knudsen et al., 1979], 

BKGI. BKGI is the current to the RPA collector measured just before the beginning of an ion sweep with 
the retarding potential set at approximately +37V relative to plasma potential. 

VPI is the value of the spacecraft potential relative to plasma potential.during an ion sweep. The 
value of VPI is derived from the least-squares fit when two or more ion masses are discemable in the 
sweep data. Otherwise, VPI is specified with a value interpolated or extrapolated from VPTE or VPI from 
nearby sweeps. 

TOTE: TOTE is the total electron density derived from the thermal electron mode saturation current FIE. 
The formula used for this present NSSDC submission, in FORTRAN language, is: 

TOTE = 6.I5E9*MAX(0, -3.5E-9 -F1E)**0.847 

We consider this measure of the total electron density to be approximate and valid only while the PV 
spacecraft is within the ionosphere. TOTE, TE, and VPTE are stored in the UADS file only if the value of 
TOTE exceeds 3E+03 cm-3. 

TE: TE is the thermal electron temperature derived using equation (1) in Knudsen et al. [1980], When 
the spacecraft is positive relative to the plasma potential, a condition existing with the spacecraft in the sun 
and in a low density plasma, the value of TE is representative of the secondary electrons trapped in the 
positive spacecraft potential well. 

FIE. FIE is the saturation electron current measured at the beginning of a thermal electron mode sweep. 
The front (retarding) grids are at a potential of +6.8 V relative to the spacecraft ground. 

BKGE: BKGE is the current measured by the RPA electrometer at the beginning of the thermal electron 
mode. The front (retarding) grids of the RPA are held at a potential of-4.6 V during the measurement. 


VPTE. VPTE is the spacecraft potential relative to the ambient plasma potential. It is derived from the 
thermal electron sweep data as described BY Knudsen et al. [1980], When the spacecraft is in the solar 
wind and exposed to the sun, its potential is typically a few volts positive with respect to the solar wind 
plasma potential. VPTE loses its meaning in this situation. 

N1: N1 is the density of the low temperature Maxwellian electron distribution used to fit the suprathermal 
electron I-V curve and is an approximation to the themal electron density. N1 is stored in the UADS file 
only when the analysis program thinks the spacecraft is within the ionosphere and assumes the vale of 
VPPE is -I.5V[Knudsen et al., 1985]. 

TI: T1 is the temperature of the low temperature Maxwellian electron distribution. 

N2: N2 is the density of the high temperature Maxwellian electron distribution used to fit the suprathermal 
electron I-V curve [IGiudsen et al., 1985], It is the density of photoelectrons, primarily, when the 
spacecraft is within the ionosphere, and is the density of solar wind electrons, primarily, when the 
spacecraft is above the ionosphere. The spacecraft is assumed to be within the ionosphere when VPPE has 
the value, - 1.5V. 

T2: T2 is the temperature of the high temperature Maxwellian electron distribution. 
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F1P FIP is the electron current measured by the RPA when operating in the suprathermal electron mode 
with zero retarding potential on the retarding grid. 

BKGP: BKGP is the electron current to the RPA when operating in the suprathermal electron mode with 
the retarding potential on the retarding grid equal to -58V. 

VPPE: VPPE is the spacecraft potential relative to the ambient plasma potential when the RPA is sweeping 
in its suprathermal electron mode. When the spacecraft is in the solar wind and not in the Venus umbra, 
the spacecraft is positive, and the potential is derived from the suprathermal electron I-V curve. When the 
spacecraft is within the ionosphere or in the Venus umbra, the potential is assumed to be -1.5V. 

The quantities TOTI, H+, 0+, 02+, C02+, TI, Nl, Tl, N2, T2 and the quantity 
sqrt(Vx**2 + Vy**2 + Vz**2) are derived by least-squares fitting a strongly non-linear numerical 
algorithm to an I-V curve. It is necessary in performing such a fit to supply an initial estimate of the 
quantities that are to be derived. If the estimates are not sufficiently close to the true least-squares values, 
the algorithm may yield a grossly erroneous value by converging to a relative minimum of the variance and 
not to the absolute minimum. Also, it may not converge at all. Although some such erroneous values have 
been eliminated from our basic tables tapes by checking for the magnitude of the variance, some erroneous 
values are known to be present. Such values can be way outside the nominal uncertainty quoted in 
Table 1. 

The algorithm that scans the data in an ion I-V curve and computes the initial estimates of the ion 
quantities must also make a decision as to what ion mass is represented by a peak in delta I (DI) [Knudsen 
et al., 1979; 1980]. The voltage at which the DI peak occurs for a given mass may be substantially smaller 
or larger than the nominal value because the Venus ionosphere is moving relative the planet with a velocity 
that can approach that of the spacecraft. Consequently, some peaks in DI have been assigned the wrong 
mass. The result is not only an erroneous concentration for that mass but also an erroneous ion velocity 
and total ion density. It is possible to recognize the incorrect assignments when comparing several I-V 
curves which are adjacent to each other in time, but the analysis algorithms are not this sophisticated. A 
few errors in ion quantities are present in the NSSDC files resulting from this difficulty. 

The analysis of a suprathermal electron I-V curve is similarly difficult. The interpretation of the electron 
distributions contributing to the I-V curve depends on the potential of the spacecraft relative to the ambient 
plasma which, in turn, depends on the location of the spacecraft and the properties of the ambient plasma. 
The spacecraft is negative in the dense ionospheric plasma. It is positive in the low density solar wind 
plasma provided the spacecraft is not in the umbra of the planet. An additional complication arises in that 
the sign of the current to the electrometer occasionally changes from negative to positive during a sweep. 
This can occur because the background current, with maximum retarding potential applied to the retarding 
grids, is compensated close to the noise level of the electrometer just before the sweep begins [Knudsen et 
al., 1979]. If the background current that has been compensated is significant relative to the saturation 
current and changes in the right direction during the ensuing sweep, the total current will go through zero 
and the sign change. The background current can change because the orientation of the rotating spacecraft 
relative to the sun changes, because a purely temporal change occurs or because the location of the 
spacecraft changes. Switching of the current from one sign to the other with the electrometer in its most 
sensitive mode produces a noise spike in the electrometer that is digitized and becomes part of the I-V 
curve. Writing an algorithm that recognizes the noise spike and the change in current sign is difficult 
because the sign of only the saturation current and background current of a sweep has been retained in the 
I-V data for reasons of minimizing the telemetry requirements of the RPA. A trained observer, looking at 
an I-V plot, can rather quickly recognize in most cases when this condition has occurred, but it is difficult 
to write an algorithm that can recognize all the possible situations and make the necessary adjustments. 

In summary, some of the quantities contained in this submittal of RPA data to the NSSDC are erroneous 
because of bad least-squares fits to the I-V curves. These bad fits have not been detected by our current 
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algorithms for reduction of the data and have not been removed by a trained observer viewing the I-V 
curves and making an educated judgment. 
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